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ABSTRACT 

We present results from an efficient, non-traditional survey to discover damped Lya (DLA) absorption 
systems with neutral hydrogen column densities Nhi > 2 x 10^" atoms cm~^ and redshifts z < 1.65. 
In the past, identification of DLA systems at z < 1.65 has been difficult due to their rare incidence 
and the need for UV spectroscopy to detect Lya absorption at these low redshifts. Our survey relies on 
the fact that all known DLA systems have corresponding Mg II absorption. In turn, Mg II absorption 
systems have been well-studied and their incidence at redshifts 0.1 < z < 2.2 as a function of the Mg II 
rest equivalent width, Wq^"^^^, is known (Steidel & Sargent 1992). Therefore, by observing the Lya fine 
corresponding to identified low-redshift Mg II systems and determining the fraction of these that are 
damped, we have been able to infer the statistical properties of the low-redshift DLA population. In an 
earlier paper (Rao, Turnshek, & Briggs 1995), we presented initial results from an archival study with 
data from HST and lUE. Now, with new data from our HST GO program, we have more than doubled 
the sample of Mg II systems with available ultraviolet spectroscopic data. In total we have uncovered 
12 DLA lines in 87 Mg II systems with Wq^'^^^> 0.3 A. Two more DLA systems were discovered 
serendipitously in our HST spectra. At the present time the total number of confirmed DLA systems at 
redshifts z < 1.65 is 23. 

The significant results of the survey are: (1) the DLA absorbers are drawn almost exclusively from the 
population of Mg II absorbers which have Wq'^'^^^> 0.6 A. Moreover, half of all absorption systems with 
both Mg II Wo^2796 g^^^j Ye II Wg^^eoo > o.5 A are DLA systems. (2) The incidence of DLA systems per 
unit redshift, udla, decreases as a function of decreasing redshift. The low redshift data are consistent 
with the larger incidence of DLA systems seen at high redshift (Wolfe et al. 1995) and the inferred low 
incidence for DLA at z = derived from 21 cm observations of gas-rich spirals (Rao, Turnshek, & Briggs 
1995). However, the errors in our determination are large enough that it is not clear if the decrease per 
comoving volume begins to be significant at z ~ 2, or possibly does not set in until z k, 0.5. (3) On 
the other hand, the cosmological mass density of neutral gas in low-redshift DLA absorbers, ^dla, is 
observed to be comparable to that observed at high redshift. In particular, there is no observed trend 
which would indicate that ^dla at low redshift is approaching the value at z = 0, which is a factor of 
~ 4 — 6.5 lower than JIdla- (4) The low-redshift DLA absorbers exhibit a larger fraction of very high 
column density systems in comparison to determinations at both high redshift and at z = 0. In addition, 
at no redshift is the column density distribution of DLA absorbers observed to fall off in proportion to 
~ N]j] with increasing column density, a trend that is theoretically predicted for disk-like systems. We 
discuss this and other mounting evidence that DLA absorption does not arise solely in luminous disks 
but in a mixture of galaxy types. 

Although we have doubled the sample of confirmed low-redshift DLA systems, we are still confronted 
with the statistics of small numbers. As a result, the errors in the low-redshift determinations of tldla 
and ^DLA are substantial. Therefore, aside from the above evolutionary trends, we also discuss associated 
limitations caused by small number statistics and the robustness of our results. In addition, we note 
concerns due to gravitational lensing bias, reliance on the Mg II statistics, dust obscuration, and the 
sensitivity of local H I 21 cm emission surveys. 

Subject headings: cosmology: observations - quasars: absorption spectra - surveys - galaxies: 
observations - galaxy formation 



1. INTRODUCTION 

QSO absorption-line systems provide one of the most 
powerful probes of the gaseous components associated 
with cosmologically distant galaxies. Except for the small 
fraction of absorption-line systems that are ejected from 
QSOs, the absorption lines reveal the presence of interven- 



ing gaseous clouds all the way up to the redshifts of the 
most distant QSOs in the Universe. Thus, independent of 
galaxy luminosity, and provided that dust obscuration ef- 
fects are small, QSO absorption lines can probe the gas in 
galaxies up to redshift z « 5, back to a time when the Uni- 
verse was less than 10% of its present age. The absorbing 
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clouds are empirically found to have neutral hydrogen col- 
umn densities in the range 10^^ < Nm < 5 x 10^^ atoms 
cm~^, and they are thought to be associated with anything 
from primordial clouds in the intergalactic medium (pos- 
sibly the weakest "Lya forest" systems) to low-luminosity 
galaxies, interacting systems, low surface brightness (LSB) 
galaxies, and the evolved gaseous spheroidal and disk com- 
ponents of the most luminous galaxies. 

Yet, despite the wealth of information on QSO absorp- 
tion lines that has been gleaned from ground-based stud- 
ies — and now from many years of observation with the 
Hubble Space Telescope (HST) — there is still much to be 
learned about the low-redshift {z < 1.65), high column 
density absorbers referred to as damped Lya (DLA) sys- 
tems. Historically these relatively rare high column den- 
sity systems were the subject of intensive searches because, 
by comparison with the Lya absorption signature pro- 
duced by the Galaxy, they were thought to represent the 
signature of the most luminous high-redshift gas-rich disk 
galaxies (Wolfe et al. 1986, hereafter WTSC86). However, 
we now know that this is not the only explanation for DLA 
systems and, in fact, the most luminous disk galaxies may 
not even dominate the population of DLA galaxies (e.g. 
Le Brun et al. 1997, Rao & Turnshek 1998). 

The DLA systems are classically defined to have neu- 
tral hydrogen column densities Nhj > 2 x 10^° atoms 
cm~^ in statistical surveys (WTSC86). There is a practi- 
cal reason for this definition although damping wings 
in a Lya absorption-line profile with low velocity disper- 
sion (10 — 20 km s~^) set in at column densities of a few 
times 10^^ atoms cm~^, they still have relatively small rest 
equivalent widths (e.g. Wq^"^^^ « 1 A). Only Lya absorp- 
tion lines with considerably larger rest equivalent widths, 
say Wq^'^^^ > 10 A, are easily recognized in low-resolution 
spectra. In turn, low-resolution spectroscopy offers an ef- 
ficient way to survey the Lya forest and search for rare 
DLA absorption lines. Thus, a 10 A rest equivalent width 
threshold is adopted in DLA spectroscopic surveys and 
this corresponds to Nhj « 2 x 10^° atoms cm~^ in the 
limit of radiative damping. 

In the past, surveys for DLA absorption have been used 
to study the distribution of neutral hydrogen at high red- 
shift (WTSC86; Lanzetta et al. 1991; Wolfe et al. 1995, 
hereafter WLFC95). These high-redshift studies indicate 
that DLA absorbers contain the bulk of the observable 
neutral gas mass in the Universe. However, spectroscopic 
detection of these absorbers using ground-based telescopes 
can only be done at redshifts z > 1.65. The detection of 
DLA at lower redshifts, where observations of Lya require 
UV data, has been difficult due to limited HST resources. 
But for A = cosmologies and qq = 0, z < 1.65 corre- 
sponds to a lookback time of w 62% of the age of the Uni- 
verse (« 77% for qo = 0.5). Thus, the determination of the 
statistics and properties of DLA systems at these redshifts 
is crucial for understanding the formation and evolution 
of galaxies. Since the occurrence of a DLA system is rela- 
tively rare and the number of QSOs studied in the UV has, 
until recently, been small, progress in this area has been 
slow. The HST QSO Absorption Line Key Project survey 
found only one DLA system with redshift z = 1.372 to- 
wards Q0935-I-417 (Jannuzi et al. 1998) in a redshift path 
of AZ = 49. Moreover, only two DLA systems have been 
confirmed from an lUE survey (Lanzetta, Wolfe, & Turn- 



shek 1995, hereafter LWT95) which originally reported the 
discovery of 14 low-redshift DLA candidate systems. One 
of these is the Key Project system towards Q0935-I-417 and 
the other is the z = 1.010 system towards EX 0302-223 
(Pettini & Bowen 1997; Boissc et al. 1998). Four of the 
14 systems had Wg^^^^ > lOA and formed the final sample 
from which they derived the low-redshift DLA statistics; 
the z = 1.010 system towards EX 0302—223 was not one 
of these four. Thus, the lUE data and the statistical re- 
sults derived from them, which were subsequently reported 
in many studies (e.g., WLFC95, Storrie-Lombardi et al. 
1996a, Storrie-Lombardi, McMahon, & Irwin 1996b, Pei, 
Fall, & Hauser 1999), should now be considered unreliable. 

Previously, we used unbiased results that could be ex- 
tracted from HST and lUE archival spectra to place lim- 
its on the properties of the low-redshift DLA population 
(Rao, Turnshek, & Briggs 1995, hereafter RTB95). In this 
paper, we discuss our recently completed HST-FOS ob- 
serving program to discover low-redshift (z < 1.65) DLA 
systems. Coupled with our earlier archival work, we have 
now completed an unbiased study of the Lya absorption 
line in 87 Mg II absorption-line systems, 12 of which we 
find to be DLA systems which formally meet the classical 
criterion of Nhi > 2 x 10^" atoms cm^^. In addition, two 
new classical DLA systems with no existing Mg II informa- 
tion were discovered serendipitously in our HST spectra. 
Consequently, our survey has greatly increased the num- 
ber of known low-redshift DLA systems. Including the 14 
classical DLA systems which were uncovered in this sur- 
vey, at the present time there are 23 published low-redshift 
{z < 1.65) DLA systems for which Nhi > 2 x 10^° atoms 
cm^^ has been confirmed on the basis of UV spectroscopy 
of the Lya absorption-line profile or 21 cm absorption. 
Although DLA systems with somewhat lower column den- 
sities can often be detected using iJS'T spectra, we have set 
our threshold column density for detection to the classical 
value so that meaningful comparisons between the statis- 
tics of these systems at low redshift and those found in 
high-redshift surveys can be made. 

In §2, we describe the survey, some ancillary data which 
we use, and the analysis which leads to the identification 
of the DLA systems. In §3, we discuss the statistical prop- 
erties of the Mg II sample. In §4, we derive the incidence 
of low-redshift DLA systems, noLAiz), and their cosmo- 
logical mass density, iloLAiz), and compare them to de- 
terminations at high redshift. In §5, we derive the DLA 
column density distribution, /{Nhi), and consider its evo- 
lution. We find that while the incidence of DLA systems 
at low redshift shows some indication of dropping with de- 
creasing redshift to a value consistent with observations of 
the H I cross-sections of local spirals, there is no clear 
evidence that their cosmological mass density drops from 
the large values observed at high redshift to the value in- 
ferred at z = (although the error bars remain large). 
Moreover, contrary to previous claims, there is evidence 
that the DLA column density distribution flattens at low 
redshift when compared to the high-redshift DLA column 
density distribution. In §6, we discuss the limitations of 
our results as well as concerns about the possible effects of 
obscuration by dust and gravitational lensing. We discuss 
the nature of the DLA absorber population in §7 and the 
conclusions are summarized in 58. 
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2. UV SPECTROSCOPIC SURVEY FOR DLA SYSTEMS 

The details of our spectroscopic survey for DLA systems 
are presented below. First we discuss the unbiased method 
of searching for DLA which allows us to derive the statis- 
tical properties of the DLA absorbers (§2.1). We then 
describe the extent to which the survey relies on archival 
observations (§2.2) and new observations (§2.3); the role of 
21 cm observations in the exclusion and selection of objects 
is also described (§2.4). The identified DLA systems are 
presented in §2.5 and an object of special interest, which 
was a DLA candidate but is now identified as a group- 
ing of Lya absorption lines on a supercluster size scale, is 
discussed in §2.6. 

2.1. Unbiased Selection of Mgll Systems for UV 
Spectroscopic Study 

The Mg II selection method for finding DLA systems 
with z < 1.65 was discussed in RTB95. In brief, since 
all DLA systems discovered so far exhibit metal-line ab- 
sorption (Turnshek et al. 1989; Lu et al. 1993; Wolfe 
et al. 1993; Lu & Wolfe 1994), the probability of find- 
ing DLA lines is greatly increased if a search is restricted 
to QSOs whose spectra show intervening metal-line ab- 
sorption. Mg II absorption lines can be detected from the 
ground with optical spectroscopy for redshifts 0.1 < z < 
2.2, and their statistical properties are known (Steidel & 
Sargent 1992, hereafter SS92). Thus, if the fraction of 
DLA systems in an unbiased sample of Mg II systems can 
be determined, then the unbiased incidence of DLA sys- 
tems at redshifts as low as z = 0.1 can be constructed from 
the known incidence of the Mg II systems. 

Our initial Mg II sample was drawn from a list of over 
360 Mg II systems in the literature with z < 1.65 (Rao 
1994). The sample of Mg II systems used to search for 
DLA absorption was constructed after excluding (1) Mg II 
absorption lines with rest equivalent width W,^^™° < 0.3 
A, (2) systems for which the predicted position of the Lya 
line fell within a broad absorption line (BAL) trough in 
the spectrum of a BAL QSO (Turnshek 1988; Turnshek 
1997a, b), (3) previously known H I 21 cm absorbers (§2.4), 
(4) systems for which the predicted position of the Lya line 
was shortward of a known Lyman series absorption edge 
(i.e. a Lyman limit), and (5) Mg II systems within 1000 
km of the QSO emission redshift. Furthermore, in 
accordance with the criterion adopted by SS92, Mg II sys- 
tems within 1000 km s^^ of each other were considered to 
l)e a single system. This rcisultcd in a sample; of 243 Mg II 
systems which could be observed for follow-up study if un- 
limited telescope time were available. These 243 systems 
define our total possible sample. 

2.2. Update on HST and lUE Archival Data 

The archival study of RTB95 included 43 Mg II sys- 
tems which satisfied the selection criteria above (§2.1) and 
which had adequate UV spectral data (see Table 1 in 
RTB95) to search for DLA or DLA candidates. Two DLA 
lines {z = 0.859 in PKS 0454-^039 and z = 0.656 in 3C336 
(1622+259)) and four candidates were identified. With 
the availability of additional archival data, the archival 
portion of the RTB95 survey has now been updated and 
modified as explained below. 

One of the four RTB95 DLA candidates, the z = 1.0096 



system towards EX 0302—223, has now been confirmed as 
a DLA system. Two candidates, the z = 0.1602 system to- 
wards Q0151-I-045 and the z = 1.2232 system towards PG 
1247-1-268, were remeasured using HST-FOS data and they 
did not meet the minimum H I column density criterion. 
The fourth RTB95 DLA candidate, the z = 0.5761 system 
towards PG 0117-1-213, was observed with HST-FOS in 
spectropolarimetric mode (see figure 3 in Koratkar et al. 
1998 for the recalibrated spectrum). We have determined 
that this candidate system is not damped, but discuss it in 
§2.6 as an object of special interest. In addition, although 
RTB95 discounted the /[/i? spectrum of the z = 1.391 sys- 
tem in Q0957-h561A as showing Nhi > 2 x 10^° atoms 
cm~^, we have now remeasured its spectrum using higher 
resolution and signal-to-noise ratio HST-FOS data and 
find it to be classically damped. Also, we find one of the 
six new archival systems, z = 0.633 in Q1209-I-107, to be 
damped. In order to treat all of the sample data con- 
sistently, we have remeasured the H I column densities 
and uncertainties for the archival DLA systems by fitting 
Voigt profiles to the data, even when the DLA systems in 
these objects were previously studied by others (see §2.5). 
These details, including six more systems that have since 
appeared in the HST archives and which we have added 
to our observed sample, are given in Table 1. Finally, we 
note that in Table 1 of RTB95 it was erroneously reported 
that the z = 0.6330 Mg II system in Q0420-014 was part 
of our sample, however there is no available UV spectrum 
of this object in the required wavelength range. Thus, this 
system is now explicitly excluded from the observed sam- 
ple. Currently, the total number of archival Mg II systems 
included in our Lya survey is 48. 

2.3. New HST-FOS Observations 

In order to extend our original archival study (RTB95) 

and improve the statistics, we carried out a large UV spec- 
troscopic survey of low-redshift Mg II systems with HST- 
FOS in an HST Cycle 6 program. The survey was op- 
timized by selecting QSOs that: (1) had more than one 
Mg II system whose Lya line could be observed in an 
FOS spectrum but that were physically unassociated, i.e., 
the difference in absorption redshifts corresponded to more 
than 1000 km s~^ and (2) were bright enough to obtain 
adequate signal-to-noise ratios in a single orbit at the ex- 
pected position of the Lya line. The 1.0-PAIR aperture 
was used to obtained all the spectra and the FOS grating 
selection was based on the expected position of the Lya 
line. The data generally have minimum signal-to-noise ra- 
tios of S/N « 11 for the high resolution G190H/RD and 
G270H/RD gratings and S/N « 15 for the low resolution 
G160L/RD and G160L/BL gratings. The QSO limiting 
V magnitude requirements for obtaining these minimum 
S/N ratios in 27 minute exposures, which was typical for 
our sample, is shown in Figure 1 for the various gratings. 
By the end of Cycle 6, before FOS was decommissioned, 
we had used 40 orbits to obtain 38 spectra of 36 QSOs 
that included the expected location of the Lya lines of 60 
Mg II systems (actually, in one case only Ly/3 and higher 
order lines were available, but this was sufficient for our 
purposes). No data were obtained during 2 orbits due to 
technical problems with HST. 

The details of this part of the Mg II sample and associ- 
ated new HST-FOS observations are given in Table 2, with 
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Table 1 

Update to the RTB95 Archival Sample 





Coordinate 








Mg II 




A /T™ TC 

Mg 1 


T~\T A C 


bource 




Designation 






^abs 












PG 0117+213 


0117+213 


16.1 


1.491 


0.5761 


0.91 


1 


+ 





Koratkar, PID 6109 


Q 0151+045 


0151+045 


16.9 


0.404 


0.1602 


1.55 


2 


+ 




Tytlcr, PID 4396 


EX 0302-223 


0302-223 


16.4 


1.409 


1.0096 


1.16 


3 


+ 


+ 


Burbidge, PID 6224 


Uooo-t-ozi 


n 'i 'i 'i 1 Q01 
Uooo-h^5zi 


i / .0 


i .zoy 


U.yooi 


0.47 


1 






Wilis, rLD 0441 


PKS 0454+039 


0454+039 


16.5 


1.343 


0.8596 


1.53 


3 


+ 


+ 


Bergeron, PID 5351 


Q 0957+561A 


0957+561A 


17.0 


1.414 


1.3911 


2.12 


4 


+ 


+ 


Burbidge, PID 5683 


Q 1209+107 


1209+107 


17.8 


2.193 


0.3930 


1.00 


5 






Bergeron. PID 5351 










0.6295 


2.92 


5 




+ 


Bergeron, PID 5351 


PG 1247+268 


1247+268 


15.6 


2.043 


1.2232 


0.48 


1 


+ 




Burbidge, PID 5095 


PG 1248+401 


1248+401 


16.1 


1.032 


0.7729 


0.76 


1 






Bahcall, PID 5664 


PKS 1327-206 


1327-206 


17.0 


1.165 


0.8500 


2.11 


6 


+ 




Tye, PID 5654 


Q 1517+239 


1517+239 


17.3 


1.903 


0.7382 


0.30 


7 






Foltz, PID 5320 


3C 336 


1622+239 


17.5 


0.927 


0.6561 


1.29 


1 




+ 


Steidel, PID 5304 


Apparent maj 


jnitude m and 


emission 


redshift 




from Veron-Cetty & Veron (1998). 





b Reference for Zabs, W^"^"^^^, and Mg I. 



^ '+' indicates the presence of, and '— ' indicates the absence of the said absorption line. is entered where no data 

are available. 

^ The Principal Investigator of the HST archival program and the proposal identification number (PID) are listed. 



References.— (1) SS92; (2) Bergeron et al. 1988; (3) Petitjean & Bergeron 1990; (4) Caulet 1989; (5) Young, Sargent, 
& Boksenberg 1982; (6) Bergeron, D'Odorico, & Kunth 1987; (7) Sargent, Boksenberg, & Steidel 1988 
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Fig. 1. — Limiting V magnitude of QSOs in the survey set by the criterion S/N = 11 for G190H and G270H spectra, and S/N = 15 for 
G160L spectra at the redshifted wavelength of Lya. For example, the spectrum of a QSO brighter than V=19.2 observed with the G270H 
grating will have S/N > 11 at the position of a 2 = 1.4 Lya line. 
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entries taken from the Vcron-Cetty & Vcron (1998) QSO 
catalog. Specifically, cohunn 1 gives the QSO's name; col- 
umn 2 gives the QSO's 1950 coordinate designation; col- 
umn 3 gives its V magnitude; column 4 gives its emission 
redshift; column 5 gives the Mg 11 absorption redshift; and 
column 6 gives the reference from which the absorption 
redshift was obtained. Columns 7, 8, and 9 give details 
of the HST-FOS observation, i.e., grating, exposure time, 
and observation date, respectively. In cases where two or 
more systems were chosen for observation with one grat- 
ing, the grating and integration time are entered alongside 
the lowest redshift absorption system. Applicable com- 
ments are noted in column 10, with an explanation in the 
footnotes. 

Unfortunately, previously unobserved Lyman limit ab- 
sorption eliminated information on the Lya absorption line 
in 21 of the 60 new Mg II systems (see column 10 of Table 
2). In addition, the G160L/RD spectrum of S4 0248+43 
was too noisy near the redshifted {z = 0.3939) Lya wave- 
length («1700 A) to determine if the Lya absorption line 
in this particular Mg II system is damped using the ex- 
isting FOB data alone. However, we have retained this 
system in our sample since its damped nature has been 
revealed through 21 cm observations (see §2.4). 

Therefore, with the aid of archival and new observations, 
it is now possible to determine the nature of the Lya line 
in a total of 87 Mg II systems. These 87 systems define our 
current unbiased sample, which we will henceforth refer to 
as the RT sample. In §3 we discuss some of the additional 
properties of this sample (e.g. other metal lines). The dis- 
cussion in the remainder of this section primarily pertains 
to the determination of which of these 87 systems have 
associated DLA lines. 

2.4. Exclusion and Selection of DLA Systems Based on 

21 cm Data 

We excluded all previously known 21 cm absorbers from 
the RT sample because we had a priori knowledge of the 
high H I column density of the absorber. Moreover, HST 
UV spectra of the background QSOs were taken by various 
researchers for the explicit purpose of studying the damped 
line corresponding to the 21 cm absorption system. The 
four excluded systems (sec RTB95) arc: z — 0.5420 in AO 
0235+164, z = 0.4369 in 3C 196 (0809+483), z = 0.3950 
in PKS 1229-021, and z = 0.6922 in 3C 286 (1328+307). 
However, concurrent with this survey, a survey for HI 21 
cm absorption in radio-loud quasars with known Mg II sys- 
tems was being carried out by W. Lane and F. Briggs (see 
Lane et al. 1998a for some initial results). Three of the 
known Mg II systems observed during our HST-FOS sur- 
vey have now been found to be 21 cm absorbers. These are 
the z = 0.3939 system towards S4 0248+43, the z = 0.3130 
system towards PKS 1127-14, and the z = 0.2213 system 
towards 01 363 (0738+313) (Lane et al. 1998a; Lane et 
al. 2000, in preparation). The Lya line in the spectrum 
of S4 0248+43 was the only Lya line which did not have 
HST-FOS data of sufficient quality to determine if the line 
was damped. However, the availability of the H I 21 
cm data, which shows that this line must be damped (see 
§2.5.2 for further details), gives our survey a certain de- 
gree of completeness since, for all observed systems which 
were selected in an unbiased manner, we have been able 
to determine whether or not they are damped. In particu- 



lar, the originally-defined sample of 243 systems included 
Mg II systems for which no information on the H I col- 
umn density was published, either from the UV Lya line 
or for the 21 cm line. Thus, all information obtained after 
the sample was assembled is acceptable, and this does not 
bias the results. If anything, the decision to exclude all 
previously known 21 cm systems might now be viewed as 
overly conservative. It could be argued that unwarranted 
exclusion of these systems might cause us to underestimate 
the incidence and mass density of DLA systems at low red- 
shifts. In the future, if H I column density information 
becomes available for all of the Mg II systems that satisfy 
our selection criteria, all of the previously identified 21 cm 
absorbers which were excluded from the sample should be 
included at that time, because they would be considered 
free of selection bias. 

Finally, we note that the z = 0.0912 DLA system to- 
wards 01 363 (0738+313), which we detected serendip- 
itously (see §2.5.5 and Rao & Turnshek 1998), has also 
been detected in 21 cm absorption (Lane et al. 1998b). 

2.5. Mg II Systems with Identified DLA Absorption 

The QSOs with DLA detections are listed in Table 3 
along with the redshift and rest equivalent width of the 
Mg II absorption lines, the redshift and column density 
of the DLA lines and their uncertainties, and the source 
of the Lya data. We have determined that 14 DLA lines 
which meet the classical criterion for being damped are 
present in our data set. Twelve of these are found in the 
87 Mg II systems and two others were found serendipi- 
tously. Redshifts for the z = 0.3941 absorption system to- 
wards S4 0248+43, the z = 0.0912 and z = 0.2212 systems 
towards 01 363 (0738+313), and the z = 0.3127 system 
towards PKS 1127—14, which were derived from more ac- 
curate H I 21 cm data, are reported without error (Lane 
et al. 1998a, b). When deriving redshifts for the remain- 
ing DLA systems, we have not corrected for wavelength 
zero point offsets that were reported in FOS data (Rosa, 
Kerber, & Keys 1998, ISR CAL/FOS-149) since we did 
not find this necessary. For the DLA systems without 
21 cm redshifts, our derived redshifts are mostly within 
1 — 2cr of the Mg II redshifts. Thus, the current data 
are sufficient to identify any DLA lines that correspond 
to metal-line absorbers. Although Voigt profile fits that 
differed by 0.0001 in redshift are clearly distinct in fits to 
DLA lines, we have reported redshift uncertainties of 0.001 
for data taken with the G190H and G270H gratings and 
0.002 for data taken with G160L gratings, since the grating 
resolutions and signal-to-noise ratios do not permit more 
accurate wavelength determinations. 

To fit the DLA lines with Voigt profiles, all of the 
pipeline-calibrated spectra were first linearly resampled in 
wavelength to 4 pixels per resolution element in order to 
match the quarter-diode substepping procedure used dur- 
ing the observations. A continuum level was fitted to each 
sp(;ctrum interactively and the spectrum was normalized 
by this continuum. In a few cases, in medium resolution 
gratings or in cases where the H I column density was 
high in a low resolution grating, we found it necessary to 
subtract a small correction (usually less than 10% of the 
continuum) before normalizing the spectrum in order to 
force the center of the damped line to zero flux. Simi- 
lar procedures were adopted by Steidel et al. (1995) and 
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Table 2 
New HST-FOS Observations 



QSO'' 


Ooordinate 




2em^ 


Zabs 




FOS 


Exposure 


Observation 


Notes'^ 










Ms II 




l-rT'a'l'lTlP' 


Time Tsl 


Date 




Q 0002-422 


0002-422 


17.2 


2.758 


0.8366 


1 


G270H/RD 


1170 


2 June 1996 


LL 










1.5413 


2 








PHL 938 


0058+019 


17.2 


1.959 


0.6128 


3 


G190H/RD 


1590 


12 July 1996 




PKS 0119-04 


0119-046 


16.9 


1.953 


0.6577 


4 


G190H/RD 


1590 


22 July 1996 




B2 0141+33 


0141+339 


17.6 


1.450 


0.4709 


5 


G160L/RD 


1580 


24 Aug 1996 




UM 366 


0143-015 


17.7 


3.141 


1.0383 


6 


G270H/RD 


1590 


5 June 1996 












1.2853 


6 










UM 675 


0150-202 


17.4 


2.147 


0.7801 


3 


G190H/RD 


1590 


12 Sep 1996 




PKS 0229+13 


0229+131 


17.7 


2.065 


0.3723 


3 


G160L/BL 


1400 


29 Aug 1996 


LL 










0.4177 


3 








LL 


S4 0248+43 


0248+430 


17.6 


1.310 


0.3939 


5 


G160L/RD 


1660 


21 Aug 1996 


poor S/N 










0.4515 


7 










PKS 0421+019 


0421+019 


17.0 


2.055 


1.3918 


5 


G270H/RD 


1600 


20 Aug 1996 












1.6380 


5 










PKS 0424-13 


0424-131 


17.5 


2.166 


1.0345 


5 


G270H/RD 


1590 


13 July 1996 


LL 










1.4080 


5 


















1.5623 


5 










Q 0453-423 


0453-423 


17.1 


2.661 


0.7255 


8 


G190H/RD 


1700 


22 May 1996 


LL 










0.9079 


8 


G270H/RD 


2440 




LL 










1.1492 


8 








LL 










1.4592 


8 








LL 


OI 363 


0738+313 


16.1 


0.6.30 


0.2213 


9 


G160L/BL 


1520 


15 May 1996 




PKS 0823-22 


0823-223 


16.2 


>0.910 


0.9103 


10 


G270H/RD 


1610 


26 June 1996 




B2 0827+24 


0827+243 


17.3 


0.941 


0.5247 


11 


G160L/RD 


1560 


28 Oct 1996 




MG 0833+1123 


0830+115 


18.0 


2.979 


0.8032 


6 


G190H/RD 


1570 


15 May 1996 


LL 


4C 13.39 


0843+136 


17.8 


1.877 


0.6064 


12 


G160L/RD 


1540 


15 May 1996 




TB 0933+733 


0933+732 


17.3 


2.525 


1.4789 


5 


G270H/RD 


1900 


6 Dec 1996 












1.4973 


5 










PKS 0952+17 


0952+179 


17.2 


1.478 


0.2377 


5 


G160L/BL 


1430 


24 Dec 1996 




SBS 0953+549 


0953+549 


17.4 


2.579 


1.0590 


5 


G270H/RD 


1790 


29 Sep 1996 


LL 










1.0642 


5 






LL 










1.2624 


5 








LL 


MARK 132 


0958+551 


16.0 


1.760 


0.2413 


3 


G160L/BL 


1670 


11 June 1996 




TOL 1035.8-27.6 


1035-276 


19.0 


2.168 


0.8230 


13 


G160L/RD 


1510 


15 June 1996 




40 61.20 


1049+616 


16.5 


0.421 


0.2251 


9 


G160L/BL 


1660 


1 June 1996 












0.3933 


9 










PKS 1127-14 


1127-145 


16.9 


1.187 


0.3130 


14 


G160L/BL 


1450 


22 June 1996 




B2 1148+38 


1148+386 


17.0 


1.304 


0.2130 


9 


G160L/BL 


1520 


31 May 1996 












0.5533 


5 










UM 485 


1213-002 


17.0 


2.691 


1.5543 


5 


G270H/RD 


1590 


2 Dec 1996 




TON 1530 


1222+228 


16.6 


2.048 


0.6681 


3 


G190H/RD 


1600 


4 July 1996 




Q 1246-057 


1246-057 


16.7 


2.224 


0.6399 


15 


G190H/RD 


1600 


12 June 1996 


LL 


40 65.15 


1323+655 


17.5 


1.624 


1.5181 


16 


G270H/RD 


1820 


9 Nov 1996 












1.6101 


16 










PG 1329+412 


1329+412 


17.2 


1.937 


1.2820 


5 


G270H/RD 


1710 


2 June 1996 












1.6011 


5 










PKS 1354+25 


1354+258 


18.0 


2.006 


0.8585 


16 


G270H/RD 


1600 


17 Nov 1996 












0.8856 


16 


















1.4205 


16 










MO 1511+103 


1511+103 


17.7 


1.546 


0.4370 


12 


G160L/RD 


1550 


25 June 1996 


LL 


40 18.43 


1540+180 


18.0 


1.661 


0.6945 


16 


G190H/RD 


1570 


6 July 1996 


LL 










0.7294 


16 








LL 










0.7936 


16 








LL 


KP 1623.7+26.8B 


1623+269 


16.0 


2.521 


0.3292 


3 


G160L/BL 


1920 


7 July 1996 


LL 










0.8881 


3 


G270H/RD 


1630 






PG 1715+535 


1715+535 


16.5 


1.940 


0.3674 


3 


G160L/BL 


1660 


8 Nov 1996 


LL 


PKS 1821+10 


1821+107 


17.3 


1.364 


0.4738 


5 


G190H/RD 


1590 


3 June 1996 


LL 










0.5702 


5 






LL 










1.2528 


5 










30 395 


1901+319 


17.5 


0.635 


0.3901 


17 


G160L/BL 


1460 


12 Aug 1996 





^ From Vcron-Cetty & Vcron 1998. 

^ The reference for the Mg II absorption line redshift: (1) Sargent et al. 1979; (2) Lanzetta, Turnshek, & Wolfe 1987; (3) Sargent, 
Boksenberg, & Steidel 1988; (4) Sargent, Young, & Bokseiiberg 1982; (5) SS92; (6) Sargent, Steidel, & Bokseiiberg 1989; (7) Womble et 
al. 1990; (8) Carswell, Smith, & Whelan 1977; (9) Boulade et al. 1987; (10) Falomo 1990; (11) Ulrich & Owen 1977; (12) Foltz et al. 
1986; (13) Jakobsen, Perryman, & Oristiani 1988; (14) Bergeron & Boisse 1991; (15) Boksenberg et al. 1978; (16) Barthel, Tytler, & 
Thomson 1990; (17) Aldcroft, Bechtold, & Elvis 1994. 

This column indicates that the nature of the Lya line associated with the Mg II absorption system could not be determined either due 
to an intervening Lyman limit system (LL), in which case there was no continuum, or because the signal to noise ratio at the expected 
position of the Lya line was very poor. 
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Table 3 

The DLA Systems Found in the Survey 



QSO 


Coordinate 




Mg II 


Lya 


AfHz/lO^u 


^jv„,/1020 


Source'' 




Designation 


^abs 


wAiivyb (A) 

^ ' 


^abs 










UM 366 


0143-015 






1.613"^ 


0.001 


2.0 


0.2 


This Program (Rao, PID 6577) 


o4 0248+43 


0248+430 


0.3939 


1.86 


0.3941 




36*^ 


Ad 


inis Program (Rao, PID 6577) 


EX 0302-223 


0302-223 


1.0096 


1.16 


1.010 


0.001 


2.3 


0.2 


Archive (Burbidge, PID 6224) 


PKS 0454+039 


0454+039 


0.8596 


1.53 


0.859 


0.001 


4.7 


0.3 


Archive (Bergeron, PID 5351) 


OI 363 


0738+313 






0.0912^^ 




15 


2 


This Program (Rao, PID 6577) 






0.2213 


0.52 


0.2212 




7.9 


1.4 


This Program (Rao, PID 6577) 


B2 0827+24 


0827+243 


0.5247 


2.90 


0.518 


0.002 


2.0 


0.2 


This Program (Rao, PID 6577) 


TB 0933+733 


0933+732 


1.4789 


0.95 


1.478 


0.001 


42 


8 


This Program (Rao, PID 6577) 


PKS 0952+17 


0952+179 


0.2377 


0.63 


0.239 


0.002 


21.0 


2.5 


This Program (Rao, PID 6577) 


Q 0957+561A 


0957+561A 


1.3911 


2.12 


1.391 


0.001 


2.1 


0.5 


Archive (Burbidge, PID 5683) 


PKS 1127-14 


1127-145 


0.3130 


2.21 


0.3127 




51 


9 


This Program (Rao, PID 6577) 


Q 1209+107 


1209+107 


0.6295 


2.92 


0.633 


0.002 


2.0 


1.0 


Archive (Bergeron, PID 5351) 


PKS 1354+25 


1354+258 


1.4205 


0.61 


1.418 


0.001 


32 


2 


This Program (Rao, PID 6577) 


3C 336 


1622+239 


0.6561 


1.29 


0.656 


0.001 


2.3 


0.4 


Archive (Steidel, PID 5304) 



* No error is reported if Zaba is taken from 21 cm absorption-line data. 

^ The Principal Investigator of the HST program and the proposal identification (PID) number axe listed. 
° The two serendipitously discovered systems for which no Mg II information is available are also listed. 
Prom 21 cm absorption-line data and assuming Tg = 700 K (Lane et al. 2000, in preparation). 



Boisse et al. (1998) to measure the column density of DLA 
lines detected in FOB spectra. Possible reasons for the 
non-zero flux at the center of a DLA line in an FOS spec- 
trum arc scattered light within the detector and incorrect 
dark background subtractions during pipeline processing. 

A Voigt damping profile, convolved with the appropriate 
Gaussian line spread function of the FOS-grating combina- 
tion, was then fitted to the data. The Dopplor parameter 
was set to 6 = 15 km s^^ for all of the fits. Note that 
for the new HST-FOS observations we used the 1.0-PAIR 
aperture which results in the line spread function FWHM 
being 0.96 FOS diodes, i.e., 3.8 pixels. This corresponds 
to typical resolutions (FWHM) of « 2.0 A for the G270H 
grating, « 1.5 A for the G190H grating, and « 7 A for 
the G160L grating. The exact resolutions in units of A 
depend on the wavelength of the damped line. 

The theoretical profiles that best fit the observed 
damped profiles were also determined interactively. In 
many cases, Lya forest lines overlapped with the DLA 
line and care had to be taken to interactively exclude them 
from the fit. In addition, the fact that some of the data 
were fairly noisy, with typical signal-to-noise ratios of « 5 
within the damped profile, made fitting a reasonable con- 
tinuum near the damped line the most challenging task. 
Since continuum placement is also somewhat subjective, 
and this usually is the largest contribution to the uncer- 
tainties, conservative estimates for errors in Nhi were de- 
termined by shifting the normalized continuum above and 
below the most likely continuum level by an offset mea- 
sured from the la (per pixel) error array of the flux spec- 
trum. The spectrum was then renormalized to the new 
levels and Voigt profiles were again fitted to the DLA 
lines in the renormalized spectra. The resulting ranges 
in column densities derived in this way are reported as 
the formal uncertainties in Nhi- To illustrate the proce- 



dure, we show examples of Voigt profile fits with high and 
low renormalized continua to the damped profile in B2 
0827+24 (Figures 2a, b). For this spectrum, the adopted 
la continuum placement error corresponds to 11% of the 
continuum; thus, the "high" and "low" continuum levels 
were set at 1.11 and 0.89 times the best fit continuum, 
respectively. The Voigt profiles that best fit the data then 
give N^'f'' = 2.2 x lO^" atoms cm^^ and N^^f = 1.8x 10^° 
atoms cm~^. Thus the column density for this system is 
reported as Nhi = (2.0 ±0.2) x 10^° atoms cm'^ (Table 3, 
see also §2.5.7). We have found that the uncertainties in 
column densities arc generally on the order of 10%; data 
noisier than that for B2 0827+24 usually result in larger 
uncertainties while higher signal-to-noise ratio data usu- 
ally result in lower uncertainties. 

Details of the results for each of the 14 identified DLA 
systems are given below. In cases where a fit to the DLA 
profile is made, the resulting redshift is reported to four 
decimal places. However, as noted above, the accuracy 
is more realistically ±0.001 for determinations made with 
the medium-resolution gratings and ±0.002 for the low- 
resolution gratings. H I 21 cm absorption redshifts are 
adopted when possible. 

2.5.1. The UM 366 (0143-015) z = 1.613 DLA Line 

No known Mg II system was reported in the literature 
at this redshift, but we discovered a DLA line at z = 1.613 
serendipitously in our G270H spectrum. We observed this 
QSO as part of the survey because it has two Mg II 
absorption-line systems at ^; = 1.0383 and z = 1.2853 
(neither is damped). Although no Mg II information is 
available for the z — 1.613 system, Sargent, Steidel, & 
Boksenberg (1989) did detect Fe II at this redshift. Eight 
percent of the continuum had to be subtracted prior to 
fitting the profile in order to force the spectrum to have 
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zero flux at the DLA line center. Figure 3 shows the nor- 
malized spectrum in the vicinity of the DLA line and the 
best-fit Voigt profile with z = 1.6128 and Nhi = 2.0 x 10^° 
atoms cm~^ overlaid. 

2.5.2. The S4 0248+43 z = 0.3941 DLA Line 

We obtained an FOS-G160L spectrum of this object to 
study the Lya lines of Mg II systems a.t z = 0.3939 and 
z = 0.4515. The Lya line of the z = 0.4515 system is not 
damped. However, the signal-to-noise ratio in our FOS 
spectrum at the position of Lya in the z = 0.3939 sys- 
tem is very poor, and the nature of this Lya line could 
not be determined. The FOS spectrum is shown in Fig- 
ure 4. However, as mentioned in §2.4, this system was 
convincingly found to exhibit HI 21 cm absorption in 
the radio continuum of the background quasar (Lane ct 
al. 2000, in preparation). Three narrow components near 
z = 0.3941 are present. They are clustered within a 40 
km s""'^ velocity interval and have optical depths between 
15 and 20%. Assuming that the 21 cm lines are optically 
thin, their strengths suggest that the total column den- 
sity is Nhi = (5.1 ± 0.6) x lO^^r^ atoms cm'^ K"! (Lane 
et al. 2000, in preparation). For a spin temperature of 
Ts = 700 K, which is generally found to be the case (Lane 
et al. 2000, in preparation), the resulting column density 
is Nhi = (3.6 ± 0.4) x 10^^ atoms cm~^. It is unfortu- 
nate that the spin temperature had to be assumed for the 
Nhi determination. This DLA line should be reobserved 
in the UV so that data on the DLA profile can be used to 
determine Nhi independently. 

2.5.3. The EX 0302- 223 z = 1.010 DLA Line 

This system was taken from the archives. Originally, it 
was chosen as a target for observation without knowledge 
of the fact that its spectrum contains a DLA absorption 
line. It has been studied by Pettini & Bowen (1997) and 
Boisse et al. (1998). A portion of the FOS-G270H nor- 
malized spectrum, overlaid with our best-fit DLA profile, 
is shown in Figure 5. Defining the level of the continuum 
between 2450 — 2500 A was problematic in making the fit 
because of the underlying O VI/Ly/3 emission line from 
the QSO. However, the core and short-wavelength wing of 
the DLA profile fit the data well. We derive a redshift 
of z = 1.0096 and a cohmm density of Nhi = 2.3 x 10^° 
atoms cm~^, in agreement with the values derived by Pet- 
tini & Bowen (1997) and Boisse et al. (1998). 

2.5.4. The PKS 0454+039 z = 0.859 DLA Line 

This system was also taken from the archives. Again, it 
was chosen as a target for observation without knowledge 
of the fact that its spectrum contains a DLA absorption 
line. It has been studied by Steidel et al. (1995). The nor- 
malized FOS-G190H spectrum in the vicinity of the DLA 
line, with our best-fit Voigt profile overlaid, is shown in 
Figure 6. We derive Nhi = 4.7 x 10^° atoms cm~^ at 
redshift z = 0.8591 with an « 6% error in Nhi- Steidel et 
al. (1995) derived Nhi = 5.7 x 10^" atoms cm with an 
error of w 5%. Our smaller value for Nhi results from the 
exclusion of the 2268 A absorption feature (a likely Lya 
forest line) from the fit, and agrees with the value which 
was derived by Boisse et al. (1998). 



2.5.5. The 01 363 (0738+313) z = 0.0912 DLA Line 

We obtained an FOS-G160L spectrum of this object to 

study the Lya line in the z = 0.2212 Mg II system of 01 
363 (0738+313) (see §2.5.6) and found the DLA system 
at z = 0.0912 serendipitously. Since this is presently the 
lowest-redshift confirmed classical DLA system known, we 
published the result in a separate contribution (Rao & 
Turnshek 1998) and the details of this system arc given 
there. The normalized spectrum and best-fit Voigt pro- 
file with redshift z = 0.0937 and Nhi = 1-5 x lO^^ atoms 
cm~^ are shown in Figure 7. This system has also been 
found to have 21 cm absorption (Lane et al. 1998b). 

2.5.6. The 01 363 (0738+313) z = 0.2212 DLA Line 

The predicted position of the Lya line in the z = 0.2212 
system of 01 363 (0738-1-313) was observed as part of our 
survey. This system is also found to be damped. It is 
the second lowest-redshift confirmed classical DLA sys- 
tem known. The details of this system appeared in Rao 
& Turnshek (1998). The normalized spectrum and best-fit 
Voigt profile with redshift z = 0.2224 and Nhi = 7.9 x 10^° 
atoms cm~^ arc shown in Figure 8. This system has also 
been found to have 21 cm absorption (Lane et al. 1998a). 

2.5.7. The B2 0827+24 z = 0.518 DLA Line 

We obtained an FOS-G160L spectrum of this object to 
study the Lya line of the Mg II system at z = 0.5247. Fig- 
ure 9 shows the normalized spectrum in the vicinity of the 
DLA line and the best-fit Voigt profile with column density 
Nhi = 2.0 x lO^" atoms cm"^ ^^j^j redshift z = 0.5176. At 
column densities this low, it is important to properly take 
into account the convolution of the Voigt damping profile 
with the low-resolution G160L line spread function. For 
this spectrum, we did not consider the possibility of sub- 
tracting a fraction of the continuum to force the fiux at the 
center of the DLA line to zero because, based on the convo- 
lution of the damping profile with the line spread function, 
some flux at line center is expected. The discrepancy be- 
tween our redshift and the Ulrich & Owen (1977) Mg II 
redshift is somewhat puzzling. We suspect that either the 
zero point wavelength correction for our G160L spectrum 
(which we ignored) was atypically large, or there was a 
wavelength calibration problem in the original Mg II data. 

2.5.8. The TB 0933+ 733 z = 1.478 DLA Line 

According to SS92, two Mg II systems are present in 
this object at z = 1.4789 and z = 1.4973. We attempted 
to investigate Lya in both of them by obtaining an FOS- 
G270H spectrum at the expected wavelengths of Lya. For- 
tunately, the expected positions of the Ly/3 and Ly7 ab- 
sorption lines were also covered. The spectrum shows a 
partial Lyman limit at 3060 A (z « 2.36), which sets in 
just to the long-wavelength side of the expected positions 
of the two Lya lines at 3013 A (z = 1.4789) and 3036 
A (z = 1.4973). The spectrum of this object is shown 
in Figure 10, with the positions of the first three Lyman 
series lines of both systems marked. The system marked 
with subscript '1' corresponds to z = 1.4789, while the 
one marked with subscript '2' corresponds to z = 1.4973. 
Note that the strong absorption feature at « 2770 A is a 
blend of Lya lines and is not damped. 
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Based on the appearance of the spectrum near Lya at 
3000—3025 A, at least one of the Lya hnes m the Mg II sys- 
tems at z = 1.4789 and z = 1.4973 appears to be damped, 
although the interpretation is confused by the Lyman limit 
system near z « 2.36. The determination of which of these 
systems is damped is better resolved with the higher order 
Lyman lines, particularly Ly/3. A normalized spectrum in 
the vicinity of the damped Ly/3 line, and its best-fit Voigt 
profile with z = 1.4784 and Nhi = 4.2 x 10^^ atoms cm~^ 
overlaid, is shown in Figure 11. We estimate the uncer- 
tainty in the derived column density to be « 20%; the 
high uncertainty reflects the low signal-to-noise ratio near 
the Ly/3 line. The same redshift and column density are 
consistent with the absorption line at the position of Ly7. 
The Ly/3 line corresponding to the z = 1.4973 system is 
not damped. 

Another interesting feature in this spectrum is the broad 
absorption depression between w 2520 A and « 2600 
A. However, this object is not a BAL QSO (Afanasev 
ct al. 1990; SS92). The absorption complex includes 
both Ly/3 lines at z = 1.4784 and z = 1.4973. Assum- 
ing that the entire absorption complex is due to Ly/3, 
then these wavelengths correspond to the redshift range 
1.457 < z < 1.535. If interpreted as a Hubble flow, 
this corresponds to « ^Thg^ Mpc along the line-of-sight 
{qo = 0.5). The extent of the absorption is similar to 
what is seen in the spectra of the QSOs TOLOLO 1037-27 
and 1038-27 (Jakobsen et al. 1986; Sargent & Steidel 1987; 
Jakobsen & Ferryman 1992), which is speculated to bo ab- 
sorbing material distributed on a supercluster size scale. 
However, high-ionization metal-line absorption (e.g. C IV) 
is present in the TOLOLO pair, but this does not appear 
to be the case here. Unfortunately, given the proximity 
of the DLA line {z = 1.478) and the Lyman limit edge at 
z » 2.36, a corresponding broad feature at the predicted 
wavelengths of Lya between 1.457 < z < 1.535 cannot be 
searched for to resolve the situation. On the other hand, 
if the 2520 — 2600 A feature is primarily due to Lya , and 
the two Ly/3 lines just happen to lie within this complex 
of absorption lines, then the feature would extend between 
1.073 < z < 1.139 and it would have an extent of « 47/ig5^ 
Mpc {qo = 0.5). 

2.5.9. The PKS 0952+17 z = 0.239 DLA Line 

We obtained an FOS-G160L spectrum near the expected 
position of Lya in the FKS 0952+17 z = 0.2377 Mg II sys- 
tem and find it to be damped. The normalized spectrum 
and our best-fit Voigt profile are shown in Figure 12. The 
fitting procedure yields Nhi = 2.1 x 10^^ atoms cm~^ at 
z = 0.2394. Eleven percent of the continuum level had to 
be subtracted from the spc;ctruni prior to the fit in order 
to force the line center flux to zero. 

There are two other absorption lines (most likely Lya ) 
on each side of the DLA line, which is suggestive of clus- 
tering on a wide velocity scale. If the Mg II absorber lies 
within a structure that includes these systems, this group 
would be an example of absorbing material on a size-scale 
which exceeds a supercluster. Interpreting the velocity in- 
terval as part of the Hubble flow would suggest that it 
extends « 92/ig5^ Mpc along the line-of-sight {qo = 0.5). 
Given the signal-to-noise ratio of the current data and the 
possibility of a clustering of lines, we believe that it is im- 
portant to reobserve this Lya line with improved signal-to- 



noisc ratio and resolution. For example, if improved data 
show that this system is simply a unique blend of many 
weak lines which resembles a DLA line (we consider this 
possibility unlikely), it would eliminate the third lowest 
redshift DLA system from our statistical sample. More- 
over, it has a moderately large column density, and thus 
plays an important role in our determination of ^Idla at 
low redshift (§4.3). Note that this absorber is not detected 
in absorption at 21 cm (Lane et al. 2000, in preparation); 
however, it would not be unique in this respect since other 
DLA lines in radio-loud quasars have gone undetected in 
21 cm absorption (e.g. Briggs & Wolfe 1983). 

2.5.10. The Q0957+561A z = 1.391 DLA Line 

In RTB95 we used lUE archival spectra to discount the 
z = 1.3911 Mg II system in this two-component gravi- 
tational lens as having a classical DLA line. Turnshek 
& Bohlin (1993) measured its column density using a co- 
addition of many low-resolution lUE spectra and obtained 
Nhi = (7 ± 2.3) x 10^^ atoms cm"^ However, with the 
higher-resolution if5'r-FOS-G270H spectra now available 
for each component, we have remeasured the H I column 
density of this system by fitting a Voigt profile to the DLA 
lines in each component of the lens. The system in compo- 
nent B has Nhi = 7 x 10^^ atoms cm~^ and thus, would 
not enter into our sample. However, we derive a column 
density of Nhi = 2.1 x 10^*^ atoms cm~^ for the system in 
component A, so it does enter into our sample, providing 
gravitational lensing does not produce a bias (see below). 
(Note that if the system in component A and component 
B had both met the minimum column density criterion, 
this would only have been counted as one system.) In per- 
forming the fits, we had to account for the Lya emission 
line from the QSO since it overlaps the long-wavelength 
wing of the DLA line. Part of the normalized spectrum of 
component A is shown in Figure 13 along with the best- 
fit DLA profile. Prior to normalizing the spectrum fit, 
we had to subtract 2.4% of the continuum level from the 
spectrum to force the center of the DLA to have zero flux. 
The DLA line and the emission lines, Lya and N V, were 
fitted simultaneously. The emission lines were fitted with 
the sum of four Gaussians added on to a continuum level 
of unity; the sum was then attenuated by a Voigt profile 
at the position of the DLA line. Each Gaussian profile 
was specified using three parameters: wavelength, ampli- 
tude, and FWHM. The number of parameters in the fit 
is large, and the best-fit parameters for the emission line 
profiles are not unique. Therefore, we ran a number of 
fits to span a reasonable range of emission-line parameters 
and determined the best-fit Voigt profile for each case after 
convolving the spectrum with a Gaussian whose FWHM 
is given by the resolution of the G270H grating. As shown 
in Figure 14, a column density Nhi = 1.6 x 10^" atoms 
cm~^ was clearly too low while Nhi = 2.6 x 10^° atoms 
cm~^ was too high, with Nhi = 2.1 x 10^" atoms cm~^ 
being the best fit. Thus, in Table 3 we report the column 
density as Nhi = (2.1±0.5) x 10^° atoms cm~^ at absorp- 
tion redshift z = 1.3907. Small variations in the emission 
line parameters had little affect on the column density de- 
termination. In fact, inclusion of the fourth Gaussian at 
2951 A, which was necessary to obtain a better fit to the 
emission line profile, had no affect on the column density 
determination because its relative amplitude was low and 
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its wavelength was far removed from the DLA Hne. Our 
results are in agreement with the measurements of Zuo et 
al. (1997) who derive Nhi = 1.9 ±0.3 x lO^" atoms cm'^ 
for component A. 

One should be careful when selecting a gravitationally 
lensed QSO for inclusion in a DLA absorption-line sur- 
vey which makes use of a magnitude limited sample of 
QSOs. This is because magnification of gravitationally 
lensed QSOs might introduce them into samples that they 
would not have otherwise entered. We discuss this fur- 
ther in §6.2.2 (see also Smette, Claeskens, & Surdej 1997). 
Our survey is susceptible to this bias for two reasons: (1) 
the archival part of the survey included data taken with 
the lUE which could observe only the brightest QSOs and 
(2) our HST-FOS observations were chosen to maximize 
the signal-to-noise ratio in one HS T oihit (see Figure 1). 
Thus, we must exclude any gravitationally lensed QSO if 
it is determined that it would not have been observed with 
either of these instruments had it not been lensed. Lensing 
models show that component A (V" = 17) has been magni- 
fied by a factor of 5.2 (Chae 1999), i.e., the unlensed QSO 
has a magnitude of 18.8. However, since Figure 1 shows 
that a QSO of magnitude 18.8 with an absorption system 
at z = 1.3911 would still have been included as a target 
for observation with HST, we have retained this system in 
our sample. 

2.5.11. The PKS 1127~U z = 0.3127 DLA Line 

We obtained an FOS-G160L spectrum of this quasar 
to study the Lya line of the z = 0.3130 Mg II sys- 
tem and found that the Lya line is damped. Figure 15 
shows the normalized spectrum in the vicinity of the DLA 
line and the best-fit Voigt profile with column density 
Nhi = 5.1 x 10^^ atoms cm'^ and redshift z = 0.3127. 
Four percent of the continuum level was subtracted from 
the spectrum before the fit to force the DLA line center to 
have zero fiux. This system is a 21 cm absorber (Lane et 
al. 1998a). 

2.5.12. The Q1209+107 z = 0.633 DLA Line 

This system was taken from the archives. It was one 
of the six new Mg II systems added to the archives after 
the initial study of RTB95. Originally, it was chosen as a 
target for observation without knowledge of the fact that 
its spectrum contains a DLA absorption line. A portion of 
the normalized FOS-G160L spectrum of this object, with 
the best-fit Voigt profile overlaid, is shown in Figure 16. 
The poor signal-to-noise ratio (w 2) near the position of 
the DLA Hne (« 1975 A) leads to a rather large (50%) un- 
certainty in the measured column density (Table 3). We 
derive Nhi = 2.0 x 10^° atoms cm"^ at z = 0.6325. This 
result is consistent with that of Boissc ct al. (1998). 

2.5.13. The PKS 135^+25 z = 1.418 DLA Line 

We obtained an FOS-G270H spectrum of this quasar 
to study the Lyo line; of the z = 1.4205 Mg II sys- 
tem and found that the Lya line is damped. Figure 17 
shows the normalized spectrum in the vicinity of the DLA 
line and the best-fit Voigt profile with column density 
Nhi = 3.2 x 10^^ atoms cm"^ redshift z = 1.4179. 
Four percent of the continuum level was subtracted from 
the spectrum before the fit to force the DLA line center to 
have zero flux. 



2.5.14. The 3C 336 (1622+239) z = 0.656 DLA Line 

This system was taken from the archives. Steidel et 
al. (1997) studied the absorption-line systems toward this 
quasar in considerable detail. However, originally it was 
chosen as a target for observation without knowledge of 
the fact that its spectrum contains a DLA absorption-line. 
A portion of the normalized FOS-G190H spectrum of this 
object and our best-fit Voigt profile is shown in Figure 18. 
We derive z = 0.6555 and Nhi ~ 2.3 x 10^" atoms cm~^ 
for the DLA system. This is consistent with the result 
found by Steidel et al. (1997) to within our quoted error 
(Table 3). 

2.6. The PG 0117+213 z = 0.576 System 

As noted in §2.2, the recalibrated HST-FOS spectropo- 
larimetric data on PG 0117+213 has been added to the 
archival portion of our sample since it covers the predicted 
wavelength of Lya in the z = 0.5761 Mg II system. The 
data were kindly provided to us by A. Koratkar. The 
spectrum is shown in Figure 19. We find that the strong 
absorption feature at « 1920 A, which appears as a can- 
didate DLA line in figure 3 of Koratkar et al. (1998), is 
more convincingly fitted with a blend of 9 Gaussian com- 
ponents. A Voigt damping profile with column density 
Nhi = 3.0 x 10^" atoms cm^^ is overplotted to illustrate 
that a single damping profile is a poor fit to this complex 
of absorption. Therefore, we have discounted it as a DLA 
line. In this system, the Lya absorption line which cor- 
responds to the published Mg II absorption redshift lies 
at 1917 A. It is the strongest line in the complex blend of 
absorption lines. This blended feature is evidently another 
example of correlated structure on a large (supercluster) 
scale. Assuming a cosmological Hubble flow, the compo- 
nents at the two extremes of the blend would be separated 
by « 34/ig5^ Mpc (go = 0.5) along the line-of-sight. 
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Fig. 2. — This figure illustrates the technique used to estimate uncertainties in Nhi- The largest source of uncertainty in measuring Nhi 
comes from continuum placement in spectra with low signal-to-noise ratios. Best fit profiles to the damped line in the spcctnun of B2 0827-1-24 
are shown here, but with continua set at levels equal to la above (panel a) and below (panel b) the most likely continiuiin level near the 
position of the damped line (Figure 9 shows the fit with the most likely continuum level). For this spectrum, la = II of the most likely 
continuum, and the column densities derived from the fits are 

j^high ^ 2.2 X 1020 atoms cm ^ from panel (a) and Ar|^]" = 1.8 x lO^" atoms 

cm~2 from panel (&). 
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Fig. 3.— Voigt profile fit to the DLA line in the normalized spectrum of UM 366 (0143-015) with z = 1.6128 and Nhi = 2.0 x 10^° atoms 
cm~2. The la error spectrum is shown as the dotted line. 
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Fig. 4. — A section of the unusable FOS-G160L spectrum of S4 0248+43 is shown with the position of the Lya line corresponding to the 
^abs = 0.3939 Mg II system marljed. The 1<t error spectrum is shown as the dotted line. 




2400 



2450 
Wavelength [A] 



2500 



Fig. 5. — Voigt profile fit to the DLA fine in the normahzed spectrum of EX 0302-223 with z = 1.0096 and Nhi = 2.3 x 10^° atoms cm'^. 
The la error spectrum is shown as the dotted line. 
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Fig. 6.— Voigt profile fit to the DLA line in the normalized spectrum of PKS 0454+039 with z = 0.8591 and Nhi = 4.7 x 10^" atoms 
cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 7.— Voigt profile fit to the DLA line in the normaUzed spectrum of OI 363 (0738+313) with z = 0.0937 and Nhi = 1.5 x lO^^ atoms 
cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 8.— Voigt profile fit to the DLA line in the normalized spectrum of OI 363 (0738+313) with z = 0.2224 and Nhi = 7.9 x 10^" atoms 
cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 9.— Voigt profile fit to the DLA line in the normalized spectrum of B2 0827+24 with z = 0.5176 and Nhi = 2.0 x lO^o atoms cm-^. 
The la error spectrum is shown as the dotted line. 
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Fig. 10. — The HST-FOS G270H spectrum of TB 0933+733. The la error spectrum is shown as the dotted line. Lyman hmit systems are 
seen at 2300 A (z = 1.522) and at 3060 A {z = 2.255). The positions of the Lyman series lines of the two Mg II absorption systems at (1) 
z = 1.4787 and (2) z = 1.4973 are marked. The absorption hne at 2770 A is a blend. Absorption in the form of a broad trough extending 
between 2520 A and 2600 A can also be seen; it suggests correlated structures on supercluster size scales (see text). 




Fig. 11.— Voigt profile fit to the damped Ly^ line in the normalized spectrum of TB 0933+733 with z = 1.4784 and Nhi = 4.2 x lO^i 
atoms cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 12.— Voigt profile fit to the DLA line in the normaUzed spectrum of PKS 0952+17 with z = 0.2394 and Nhi = 2.1 x lO^^ atoms 
cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 13. — Composite profile fit to the emission lines Lya and N V, and the DLA absorption line in the spectrum of Q0957+561A. The 
emission lines are best fit by a sum of 4 Gaussians added on to a continuum level of unity. The best fit Voigt profile has z = 1.3907 and 
Nhi = 2.1 X 10^° atoms cm . The la error spectrum is shown as the dotted line. 
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Fig. 14. — The DLA absorption line in the spectrum of Q0957+561A with the best fit composite profile shown as the solid line 



[Nhi = 2.1 X lO^O atoms cm"^), 

is clearly too low and Nhi = 2.6 X 10"^" atoms cm""^ (lower dashed line) is clearly too high. We report the column density of this system as 
Nhi = (2.1 ± 0.5) X 10^0 atoms cm . The la error spectrum is shown as the dotted line. 



The dashed lines show that a column density of Njjj = 1.6 X 10 atoms cm (upper dashed line) 
10^" 




Fig. 15. — Voigt profile fit to the DLA fine in the normalized spectrum of PKS 1127—14 with z = 0.3127 and Nhj = 5.1 x 10^^ atoms 
cm~^. The la error spectrum is shown as the dotted line. 
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Fig. 16.— Voigt profile fit to the DLA line in the normalized spectrum of Q1209+107 with z = 0.6325 and Nhi = 2.0 x lO^o atoms cm-2. 
The 1<T error spectrum is shown as the dotted line. 
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Fk;. 18. Voigt profile fit to the DLA fine in the normahzed spectrum of 3C 336 (1622+239) with z ■■ 
■xa~^ . The 1<t error spectrum is shown as the dotted fine. 
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Fig. 19. — The strong absorption feature in the H ST-FOS G190H spectrum of PG 0117+213 is fit with a blend of 9 Gaussian components 
(dotted line). A single Voigt damping profile with Nhi = 3.0 X 10^'^ atoms cm~^ is overplotted as the dashed line and is a poor fit to the 
data. The deepest component is at the redshift of the Mg II system and has rest equivalent width Wq^^^^ = 2.7 A. The la error spectrum 
is shown as the dot-dashed line. 



19 



3. STATISTICAL PROPERTIES OF THE Mg II SAMPLE 

Table 4 lists the results on the detection of DLA in the 

87 Mg II systems which make up our low-rcdshift sample. 
We call this sample the RT sample to distinguish it from 
the SS92 Mg II sample. Column 1 of Table 4 gives the 
QSO's 1950 coordinate designation; column 2 gives the 
absorption redshift; columns 3 — 6 give the available rest 
equivalent widths (or upper limits) of the Mg II A2796, 
Mg II A2803, Mg I A2852, and Fe II A2600 absorption 
lines, respectively; column 7 gives the references for the 
absorption-line data; and column 8 gives the column den- 
sity of the classical DLA absorption line if present in the 
system. 

In RTB95, we discussed the rationale behind using a 
sample of Mg II systems with a boot-strapping method to 
determine the statistical properties of low-redshift DLA 
systems. Inherent in this method is the reliance on the 
empirically determined statistical properties of Mg II sys- 
tems by SS92 and the assumption that the Mg II samples 
of SS92 and RT represent the same parent population of 
absorbers for Wi^'^'^^^ > 0.3 A. Therefore, to justify the 
results which we present in §4, we discuss the statistical 
properties of the SS92 and RT samples and show that they 
are equivalent (§3.1). We also discuss the properties of 
classical DLA absorbers in the RT sample and show that 
they are primarily drawn from the population of Mg II ab- 
sorbers with Wo^2796 > 0.6 A (also §3.1). We consider the 
possible presence of trends between the existence of clas- 
sical DLA absorption and: (1) the Mg II absorption-line 
doublet ratio (§3.1.1), (2) Fe II A2600 absorption (§3.1.2), 
(3) and Mg I A2852 absorption (§3.1.3). Finally, we derive 
the effective redshift path covered in our survey (§3.2). 

3.1. The Mg II A2796 Rest Equivalent Width Distribution 

Figure 20 shows the Mg II A2796 rest equivalent width 
(Vl/o^2796) distributions for the SS92 sample (A^ss92), the 
RT sample {Nut), and the DLA subsample of the RT sam- 
ple (Nrt-dla)- The SS92 Mg II data have been scaled 
to the RT Mg II data. It is clear that the SS92 and RT 
Wo^2796 distributions are very similar; this is demonstrated 
more rigorously below. It is also striking that while 43 of 
the 87 RT systems have 0.3 < Wo^^rge ^ q q ^^i^ 
of the 12 DLA systems falls in this interval. The remain- 
ing 11 DLA systems have Wp^^Tge > q q ^ Therefore, the 
probability of finding a DLA system in a Mg II sample 
with Wo^2796 > 0.3 A is 12/87 = 13.8%, but the proba- 
bility becomes 11/44 = 25.0% if the Wo^^rge threshold is 
increased to VKo^2796> o.6 A. 

The normalized cumulative Wq'^''^^^ distribution func- 
tions for the three above samples (SS92, RT, and RT-DLA) 
are shown in Figures 21 and 22 for the Wq^^™^ > 0.3 A and 
^A27g6 > 0.6 A thresholds, respectively. The maximum 
deviations between these distributions have been used in 
a Kolmogorov-Smirnov (KS) test to determine the prob- 
ability that any two of the samples are drawn from the 
same parent population. The results of the KS tests are 
given in Table 5. The SS92 and RT samples are com- 
pletely consistent with being drawn from the same parent 
population of absorbers for Wq'^'^^^> 0.3 A; and the sam- 
ples are nearly identical for W(^^^''*'> 0.6 A. This validates 
our assumption that the statistical properties of the SS92 

^However, this may not be true for the population of weak Mg II 



Mg II sample arc; applicable to the RT Mg II sample. On 
the other hand, the KS statistic shows that there is only 
a 2% probability that the RT-DLA subsample is drawn 
from the same parent population as the Wq"^"^^^ > 0.3 A 
RT sample, however the probability increases to 36% for 
the W^'^''^^ > 0.6 A RT sample. Aside from this thresh- 
old correlation (the tendency for Nhi > 2 x 10^° atoms 
cm~^ systems to be present when Wq'^'^^^> 0.6 A), Figure 
23 shows that there is no further trend between Nhi and 
T-Vq^^^^^ above these thresholds. Nevertheless, wc suspect 
that this threshold correlation is revealing something fun- 
damental about the absorbers that give rise to classical 
DLA absorption-line systems. 

Ionization equilibrium studies have shown that QSO 
absorbers must be optically thick at the Lyman limit 
{Nhi > 3 X 10^'' atoms cm~^) for appreciable amounts 
of Mg+ to be present (Bergeron & Stasiriska 1986). This 
is because Mg+ will only survive the presence of a strong 
background radiation field if it is embedded in an opti- 
cally thick H I cloud which will shield it from the ion- 
izing radiation. This prediction has been verified since it 
has also been empirically shown that Mg II and Lyman 
limit systems arise in the same population of absorbers 
(SS92 and references therein). Thus, it is understood that 
Nhi > 3 X 10"'^'' atoms cm~^ is generally necessary for the 
detection of a Mg II system"'. Our above empirical analysis 
suggests that, in addition, Wi^^'^^^> 0.6 A is almost always 
necessary (but not sufficient) to have Nhi > 2x 10^" atoms 
cm~^. Why should this be the case? Before we discuss this 
(§3.1.4), we wish to note the empirical results from several 
other comparisons which we have considered in an attempt 
to understand the nature of the DLA population. 

3.1.1. The Mgll Doublet Ratio and DLA Absorption 

Figure 24 is a plot of the Mg II doublet ratio, 
M/>27g6/|^A2803^ vcrsus Wo^27g6 for the RT sample. The 

DLA systems from the RT sample are marked with open 
squares. We have also plotted measurements from the 

literature for the four 21 cm absorbers that we excluded 
from our unbiased Mg II sample (§2.4). Although we have 
not included these previously known 21 cm absorbers in 
the statistical analysis of low-redshift DLA systems, their 
properties may be relevant for assessing the nature of the 
low-rcdshift DLA population. The horizontal dashed line 
in the figure has W^^''^^ /Wff^°^ = 1.5; all of the DLA 
systems have Mg II doublet ratios that fall below this line 
indicating that their Mg II lines lie somewhat close to the 
saturated part of the curve of growth. However, the Mg II 
doublet ratio exceeds 1.5 about 30% of the time in non- 
DLA systems. 

3.1.2. Mgll A2796 versus Fell A2600 and DLA 
Absorption 

Bergeron & Stasinska (1986) used photoionization mod- 
els to show that Mg II systems with strong Fe II absorp- 
tion have H I column densities in excess of a few times 
lO'^^ atoms cm~^. Le Brun et al. (1997), in their search 
for DLA galaxies, used the criterion Wo^27g6/|^A260o ^ ^ 
as a predictor of DLA in an absorption-line system. 



absorbers with Wq^2796< q.sA (Churchill et al. 1999a). 
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Fig. 20.— The Mg II A2796 rest equivalent width (11/(5^2796) distributions for the SS92 sample (A'ss92). the RT sample (Nrt), and the 
DLA subsample of the RT sample (AfflT— dla)- The SS92 Mg II data have been scaled to the RT Mg II data. The figure shows that the 
SS92 and RT distributions are very similar. Also note that while only 1 of the 43 RT systems with 0.3 < Wq^"^^^ < 0.6 A is a DLA system, 
11 of the 44 RT systems with Wg^^Tse > q q ^ gj-g j^LA systems. 
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Fig. 21. The normalized cumulative Mg II rest equivalent width {W^^^'^^^) distributions for the SS92, RT, and RT-DLA samples with 
(^A279(> y Q 3 ^ 'j'j^g figure shows that the RT and SS92 distributions are nearly identical. The maximum deviations between any two of 
these distributions are used in a K-S test to determine the probability that they are drawn from the same parent population. The results are 
given in Table 5. 
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Fig. 23. — A plot of the logarithm of the H I column density versus Mg II A2796 rest equivalent width for the 12 DLA systems in the RT 
sample showing that no trend exists between these two quantities. 
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Table 4 
The RT Mg II Sample'' 



QSO 


^abs 


^0^2796 (A) 


^yA2803 (A) 


^^A2852 (A) 


^A2600 (A) 


References 


iVif//1020 




Mg II 


Mg II 


Mg II 


Mgl 


Fell 




atoms cm~2 


0002-422 


1.5413 


0.48 


0.32 


<0.3 


<0.1 


1 




0002+051 


0.8516 


0.92 


0.84 


0.16 




2 




0058+019 


0.6128 


1.70 


1.51 


0.35 


1.39 


3 




0117+213 


0.5761 


0.91 


0.93 


0.23 


0.88 


2 






1.0473 


0.47 


0.26 


<0.2 


<0.1 


2 




0119-046 


0.6577 


0.30 


0.22 


<0.1 


<0.1 


4 




0141+339 


0.4709 


0.78 


0.65 


<0.3 


<0.7 


2 




0143-015 


1.0383 


0.64 


0.53 


<0.1 


<0.1 


5 






1.2853 


0.56 


0.33 


<0.1 


<0.1 


5 




0150-202 


0.7800 


0.36 


0.21 




<0.3 


3 




0151+045 


0.1602 


1.55 


1.55 


<0.5 


<0.5 


6 




0215+015 


1.3447 


1.93 


1.64 


0.32 


1.36 


7 




0248+430 


0.3939 


1.86 


1.42 


0.70 


1.03 


2 


51 




0.4515 


0.34 


0.29 


<0.5 


<0.1 


8 




0302-223 


1.0096 


1.16 


0.96 


0.18 


0.63 


9 


2.3 


0333+321 


0.9531 


0.47 


0.33 


<0.3 


<0.2 


2 




0421+019 


1.3918 


0.35 


0.31 


<0.3 


<0.1 


2 






1.6380 


0.35 


0.28 


<0.2 


<0.4 


2 




0424-131 


1.4080 


0.55 


0.35 


<0.3 


0.44 


2 






1.5623 


0.38 


0.39 


<0.2 


<0.2 


2 




0454-220 


0.4744 


1.44 


1.31 


0.35 


0.97 


10 






0.4834 


0.45 


0.27 


0.09 


0.16 


10 




0454+039 


0.8596 


1.53 


1.40 


0.37 


1.11 


2 


4.7 




1.1538 


0.57 


0.36 


<0.2 


<0.1 


2 




0710+119 


0.4629 


0.62 


0.29 


0.24 


<0.4 


11 




0735+178 


0.4240 


1.32 


1.03 


0.18 


0.87 


12 




0738+313 


0.2213 


0.52 


0.50 


<0.2 


<0.6 


13 


7.9 


0742+318 


0.1920 


0.33 


0.23 


<0.2 




13 




0823-223 


0.9103 


1.15 


0.68 


<0.4 




14 




0827+243 


0.5247 


2.90 


2.20 




1.90 


15,16 


2.0 


0843+136 


0.6064 


1.08 


0.67 


<0.1 


<0.7 


17 




0933+732 


1.4789 


0.95 


1.15 


<0.3 


0.76 


2 


42 




1.4973 


1.71 


1.98 


0.67 


1.15 


2 




0952+179 


0.2377 


0.63 


0.79 


<0.4 




2 


21 


0957+561A 


1.3911 


2.12 


1.81 


0.19 


1.67 


18,19,20 


2.1 


0958+551 


0.2413 


0.55 


0.38 


<0.2 




3 




1035-276 


0.8242 


1.08 


0.87 


<0.6 


0.52 


21 




1038+064 


0.4416 


0.66 


0.57 


<0.2 


<0.2 


2 




1040+123 


0.6591 


0.58 


0.42 


<0.2 


<0.2 


11 




1049+616 


0.2255 


0.51 


0.56 


<0.1 


<0.4 


13 






0.3937 


0.34 


0.29 




<0.1 


13 




1100-264 


1.1872 


0.51 


0.28 


<0.2 


<0.2 


9 






1.2028 


0.54 


0.43 


0.27 


<0.2 


9,22 




1115+080 


1.0431 


0.31 


0.18 


<0.2 


<0.1 


2 




1127-145 


0.3130 


2.21 


1.90 


1.14 


1.14 


23 


51 


1137+660 


0.1164 


0.50 


0.53 


<0.2 


<0.2 


24 




1148+386 


0.2130 


0.96 


0.37 


<0.3 




13 






0.5533 


0.92 


0.99 


<0.2 


<0.2 


2 




1206+459 


0.9277 


1.00 


0.79 


<0.2 


<0.1 


2 




1209+107 


0.3930 


1.00 


0.54 


<0.2 


<0.4 


25 






0.6295 


2.92 


2.05 


<0.4 


1.50 


25 


2.0 


1213-002 


1.5543 


2.09 


1.65 


<0.2 


1.11 


2 




1222+228 


0.6681 


0.43 


0.41 


<0.1 


<0.1 


3 




1229-021 


0.7571 


0.52 


0.48 


<0.1 




11 




1241+176 


0.5506 


0.48 


0.37 


<0.2 


<0.2 


2 




1246-057 


1.2015 


0.90 


0.75 


<0.4 


0.40 


2 






1.6453 


0.52 


0.58 


<0.3 


0.38 


2 
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Table 4 — Continued 



QSO 


^abs 


^A2796 (^J 


^yA2803 j'^) 


^yA2852 


pyA2600 


References Nhi/W^° 




Mg II 


Mg II 


Mg II 


Mgl 


Fell 


atoms cm~^ 


1247+267 


1.2232 


0.48 


0.38 


0.11 




2 


1248+401 


0.7729 


0.76 


0.50 


<0.4 


0.32 


2 


1254+047 


0.5193 


0.46 


0.38 


<0.2 


0.40 


2 


1317+277 


0.2887 


0.33 


0.31 


<0.2 




2 




0.6596 


0.49 


0.33 


<0.2 


0.17 


2 


1323+655 


1.5181 


0.57 


0.53 


<0.1 




26 




1.6101 


2.20 


1.85 


0.16 


0.88 


26 


1327-206 


0.8526 


2.11 


1.48 


<0.4 


0.76 


27,28 


1329+412 


1.2820 


0.49 


0.31 


<0.3 




2 




1.6011 


0.70 


0.35 


<0.2 


<0.2 


2 


1331+170 


1.3284 


1.73 


1.14 


<0.1 


0.32 


2 


1338+416 


0.6213 


0.31 


0.17 


<0.3 


<0.1 


2 


1354+195 


0.4563 


0.89 


0.82 


0.16 


0.32 


2 


1354+258 


0.8585 


1.00 


0.86 


<0.1 


<0.2 


26 




0.8856 


0.81'' 


0.57'' 


<0.2 


<0.2 


26 




1.4205 


0.61 


0.50 


0.20 


0.55 


26 32 


1421+330 


1.1725 


0.53 


0.40 


<0.1 


<0.1 


2 


1517+239 


0.7382 


0.30 


0.34 




<0.1 


3 


1622+239 


0.6561 


1.29 


<1.69 


<0.4 


1.13 


2 2.3 




0.8914 


1.55 


1.27 


<0.3 


1.32 


2 


1623+269 


0.8881 


0.93 


0.75 


0.14 


0.21 


3 


1634+706 


0.9903 


0.58 


0.42 


<0.2 


<0.1 


2 


1704+608 


0.1634 


0.59 


0.52 


<0.5 


<0.2 


13 




0.2220 


0.55 


0.33 


<0.2 




13 


1821+107 


1.2528 


0.71 


0.48 


0.09 


<0.3 


2 


1901+319 


0.3901 


0.45 


0.15 


<0.1 




11 


2128-123 


0.4299 


0.41 


0.37 


0.10 


0.27 


9 


2145+067 


0.7906 


0.52 


0.41 


<0.2 


<0.1 


2 


2223-052 


0.8472 


0.65 


0.42 


<0.2 


<0.1 


29.30 


232(1-177 


1.2010 


0.50 


0.38 


<0.1 


<().! 


1) 



^ Upper limits arc cither given by the authors or estimated by us from the pubUshed spectra. If the hne is part of a blend, then the rest 
equivalent width entered is an upper limit. The upper limit is a la upper limit if estimated by us from the published spectrum. In some 
cases, the absorption line (A2600 or A2852) was clearly visible in a spectrum but was not identified by the authors since the rest equivalent 
width of the line did not meet their detection criterion, usually 5<t. In these cases, we estimated the equivalent width of the line from the 
spectrum and tabulated the measurement as an upper limit. The trends that we have established in Figures 24, 25, and 26 are evident from 
the measurements of rest equivalent widths by the authors of the published spectra; our upper limit measurements only serve to provide a 
certain degree of completeness to our sample and, in fact, do not indicate or contribute to any trends by themselves. 

'' The rest equivalent widths quoted in reference (26) are in error. 



References. (1) Lanzctta, Turnshek, & Wolfe 1987, (2) SS92, (3) Sargent, Boksenberg, & Steidel 1988, (4) Sargent, Young, & Boksenberg 
1982, (5) Sargent, Steidel, & Boksenberg 1989, (6) Bergeron et al. 1988, (7) Bergeron & D'Odorico 1986, (8) Womble et al. 1990, (9) Petitjean 
& Bergeron 1990, (10) Tytler et al. 1987, (11) Aldcroft, Bechtold, & Elvis 1994, (12) Boksenberg, CarswcU, & Sargent 1979, (13) Boissc et al. 
1992, (14) Falomo 1990, (15) Wills 1978, (16) Ulrich & Owen 1977, (17) Foltz et al. 1986, (18) Weymann et al. 1979, (19) Caulet 1989, (20) 
Wills & Wills 1980, (21) Dinshaw & Impcy 1996, (22) Boisse & Bergeron 1985, (23) Bergeron & Boisse 1991, (24) Bahcall et al. 1993, (25) 
Young, Sargent, & Boksenberg 1982, (26) Barthel, Tytler, & Thomson 1990, (27) Kunth & Bergeron 1984, (28) Bergeron, D'Odorico, & Kunth 
1987, (29) Le Brun et al. 1993, (30) Miller & French 1978. 
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Wc have used the Mg II systems in the RT sample with 
p^A2600 measurements to further test this. Figure 25 shows 
Wq^^ versus Wq^^^^" for systems for which Fe II equiv- 
alent width measurements are available in the literature, 
or for which we were able to derive upper limits from pub- 
lished spectra (see Table 4). As in Figure 24, the DLA 
systems in the RT sample are marked with sqiiarcs and the 
21 cm absorbers which are not part of the RT sample are 
marked with circles. The horizontal and vertical dashed 
lines identify the region in the Wq'^'^^^ versus Wq'^^'^'^ plane 
for which W'o^^Tge^ q.S A and V^o^^eoo > o.5 A. Rehable 
equivalent widths are available in the literature for 24 out 
of 27 of these systems, where detections of Fe II A2600 
usually exceed the 5(t level of significance. Several results 
of interest emerge from this plot. First, there is a lin- 
ear correlation between Wq^"^^^ and Wq^^°^. If we force 
a least squares fit to the 24 data points to pass through 
the origin, we obtain Wg^^rge^ 45|;f/A2600 rpj^^ correla- 
tion between Mg II and Fe II in high-redshift, predomi- 
nantly low-excitation, absorption-line systems was noted 
by Bergeron & Stasihska (1986). They found that a strong 
Mg II system with Wq'^^^^ > 1 A was invariably accompa- 
nied by strong Fe II absorption. Second, excluding upper 
limits, 10 of the 21 systems in the RT sample (« 48%) 
with Wo^2796> 0.5 A and Wo^^eoo > o.5 A are DLA sys- 
tems. Third, the non-DLA systems in this regime all have 
H I column densities Nhi > 10^^ atoms cm^^. Fourth, 
with the exception of the known z = 0.692 21 cm ab- 
sorber towards 3C 286 which has unusually low metal-line 
equivalent widths {W^^'^^^= 0.39 A and W^Q^^eoo ^ o.22 
A, Cohen et al. 1994), all of the DLA systems lie in the 
upper-right region of the plot where Wq^'^^^> 0.5 A and 
^A2600 y 5 ^ ^j^g future, these results can be used 
to further improve the selection criteria for searching for 
DLA. In particular, there is approximately a 50% probabil- 
ity of finding a classical DLA system in metal-line systems 
with Wo^2796> 0.5 A and Wo^^eoo > o.5 A. 

3.1.3. M^rll A2796 versus Mgl A2852 and DLA 
Absorption 

Figure 26 shows W^'^''^^ versus Wg'^^^^. It can be seen 
that the DLA systems do not lie in a special region in 
the Mg II A2796 versus Mg I A2852 rest equivalent width 
plane. Evidently the presence of DLA does not guarantee 
the presence of detectable Mg I A2852 in medium resolu- 
tion spectroscopic surveys. However, although the statis- 
tics are small, all systems with Wq'^^^'^ > 0.7 A have DLA. 

3.1.4. Discussion 

The Mg II absorption lines under study in these samples 
usually have Mg+ column densities which place them on 
or near the saturated regime on a curve-of- growth plot. 
The DLA systems generally appear to exhibit somewhat 
greater saturation. In either case, VFq^^™*' is a poor indi- 
cator of Mg II column density, and a better indicator of 
the number of clouds along the sight-line and their kine- 
matic spread in velocity (e.g. Petitjean & Bergeron 1990, 
SS92, Churchill 1997). Thus, our results indicate that 
the low-redshift DLA systems predominantly arise along 
sight-lines in which the number of clouds and their spread 
in velocity is systematically larger in comparison to the 
0.3 <Wo^2796^ 0.6 A systems. The velocity spread, i.e., 



the difference between the minimum velocity and maxi- 
mum velocity components, must be > 64 km s^^, since a 
Mg II A2796 absorption line which is completely black over 
the line profile must be spread over a 64 km s~^ interval 
to produce a line with W^'^'^'^^ = 0.6 A. Empirically, since 
we find that the Mg II lines are not necessarily completely 
black, the velocity spread must typically exceed this value. 
However, observations in the HI 21 cm absorption line 
of DLA systems usually show that the components which 
give rise to the Voigt damping profile are spread over a 
much smaller velocity interval, which indicates that only 
a small fraction of the clouds contributing to Wq'^'^^^ give 
rise to the DLA line. 

Other interesting correlations which pertain to high col- 
mim density systems (but not necessarily DLA systems) 
have been identified by Churchill et al. (1999b) using high- 
resolution (6 km s~^) Keck observations of Mg II systems. 
(High-resolution data are not generally available for our 
systems.) In their figure 2b Churchill et al. (1999b) plot 
|;j/A2796 v^j-gyg cj;^2796 ^ -y^rherc uj^^'^^^ is a velocity weighted 
optical depth which is sensitive to high-velocity "outlier" 
clouds with small rest equivalent widths. The systems 
with Nhi > 10^^ cm^^ in this plot (i.e. neutral-hydrogen- 
rich systems), which include all of the systems in Fig- 
ure 25 with 1^0^2796 J, 0.5 A and Wo^^eoo > o.5 A and 
of which the DLA systems form a subset, are isolated 
and occupy the top-central region of this plane, i.e., they 
have W(^'^''^^'> 0.9 A and intermediate values of lu^^'^^''. 
Figure 2a of Churchill et al. (1999b) shows that these 
neutral-hydrogen-rich systems lie in an isolated region of 
the Mg II VFo^2796 versus C IV Wq'^^'^^ plane, corresponding 
to smaller Wq^^^^ /W^^'^^^ ratios. Lastly, figure 3b from 
Churchill et al. (1999b) shows that for the few cases in 
which the absorbing galaxies have been identified, the im- 
pact parameters are among the smallest observed. These 
are evidently positive correlations for neutral-hydrogen- 
rich systems that may also hold for classically damped 
systems. 

If we consider these pieces of information together, the 
picture that emerges suggests that the metal-line absorp- 
tion components associated with a DLA system, which 
generally correspond to components other than those 
which give rise to the DLA line itself, are systematically 
different in several ways: (1) There is a relatively large 
number of absorbing components spread over a substan- 
tial velocity interval (to produce Wq^''^^^> 0.6 A), but for 
a fixed W^^^^^ the absorbing components making up the 
associated metal-line system will be more closely spaced 
in comparison to cases where DLA is not present. This 
is suggestive of a clustering of clouds in the vicinity of 
the cloud that gives rise to the DLA line itself. (2) The 
C IV Wo^i5*VMg II W^'^'^^^ metal-hne ratio is relatively 
small and the Fe II W^^^^'^/Mg II T4^o^2796 metal-line ra- 
tio is relatively high, signifying a lower ionization phase of 
the DLA absorbing region. (3) The transverse separation 
(impact parameter) from the sight-line for any identified 
DLA absorbing galaxy is relatively small in comparison 
to systems which do not have a DLA line. This qualita- 
tive description of the kinematic and ionization structure 
of DLA systems and the associated absorbing galaxies ap- 
pears to be roughly consistent with the observations, and 
should be incorporated into any successful model describ- 
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Table 5 
K-S Test Results'" 



Samples 


Wo > 0.3 A 


Wo > 0.6 A 


RT & SS92 


0.65 


0.83 


SS92 & RT-DLA 


0.01 


0.15 


RT & RT-DLA 


0.02 


0.36 



^ Probability that the two samples are drawn from the same 
parent population. 




W, 



X2796 



(A) 



Fig. 24.— A plot of the Mg II doublet ratio, VyQ^2796/^yA2803^ versus W^^^Tge f^j. (.^g j^-p sample (see Table 4). The DLA systems from 
the RT sample are marked with open squares. The open circles represent previously known 21 cm absorbers that were excluded from our 
unbiased Mg II sample. The horizontal dashed line has w^'^'^^s /w^'^^os _ g r^^ie upward pointing arrow indicates an upper limit to the 
measured value of W,^^®"^. 
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Fig. 25. — A plot of Mg II W^ '^"^'^^' versus Fe II lVf^2600 rpjj^ j^j^j^ systems from the RT sample arc marked with open squares. The open 
cireles represent previously known 21 cm absorbers that were excluded from our unbiased Mg II sample. Left pointing arrows indicate upper 
limits to the measured value of Wq^^^^ (see Table 4). The horizontal and vertical dashed lines identify the region for which W^'^'^'-"^> 0.5 A 
and VKo^2600 > 0.5 A. All the non-DLA systems in this regime have H I column densities Njjj > 10^^ atoms cm ^. The solid line has a slope 
of 1.45. 




w, 



A2852 



(A) 



Fig. 26.— A plot of Mg II W^^^TQe versus Mg I ^0^^852 r^^ie DLA systems from the RT sample are marked with open squares. The open 
circles represent previously known 21 cm absorbers that were excluded from our unbiased Mg II sample. Left pointing arrows indicate upper 
Umits to the measured value of W^"^^^"^ (see Table 4). 
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ing the nature of DLA absorbers. 

3.2. The Effective Redshift Path of the RT Survey, /S.ZjiT 

In conventional surveys, like the SS92 Mg II survey, 
searches for absorption lines are carried out by observ- 
ing the continuum regions of an unbiased sample of QSOs 
down to some minimum equivalent width threshold over 
the total redshift path of the survey. We use the param- 
eter AZiWmin) to represent the total redshift path cov- 
ered down to absorption- line rest equivalent width, Wmin, 
at mean redshift, < z >. The AZ{Wmin) parameter is 
computed by adding all of the individual redshift inter- 
vals covered by all survey spectra with data good enough 
to detect Wmin at the A — 5a level of significance. This 
can also be expressed as the integral of the redshift path 
density, g{z,Wmin)-, over redshift, 



absorption-line systems is given by 



^Z{Wmin) 



-f 

Jo 



g{z, Wmin)dz. 



(1) 



Descriptions of these standard techniques have been given 
by WTSC86, Lanzetta, Turnshek, & Wolfe (1987), SS92, 
and others. However, the RT survey is not a conventional 
survey that samples redshifts paths along random sight 
lines. We search for DLA lines at specific redshifts along 
predetermined sight lines. Nevertheless, we do have indi- 
rect information which allows us to infer an effective total 
redshift path for Mg II A2796 covered in the RT survey, 
AZuriW^jl^^), at mean redshift, < z >rt- Note that 
this is also the redshift path for the detection of DLA 
lines. 

To determine AZuxiW};^™), we need to determine 

the effective redshift path density of the RT survey, 
9rt{z,W^^^^). In turn, 9rt{z,W^^^) can be written 
as 



(2) 



yRT(Z, yVmin ) / tj/A2796\ ' 

"■V^' min ) 

where mRT{z, W^^J^*^) is the redshift distribution of Mg II 
systems in the RT sample and n{z, W^J^^) is the number 
of Mg II systems per unit redshift as a function of redshift. 
The mRT{z, W^fJ'^) distributions for VK^^rae ^ q.S A and 
0.6 A arc shown in Figure 27. The redshift evolution of 
the number density of Mg II systems can be parameterized 

by 

n{z, W,^Z'') = no(W^™)(l + z^^^'^ (3) 

where no{W^^J^^) and ^/{W^J^^) are obtained from Table 
4 of SS92. These values are n{z = L12, W^^rae ^ q.S) = 
0.97 [which implies noiW^Z^^ = 0.3) = 0.54], n{z = 
1.17,W^2796 ^ 0.6) = 0.52 [which implies no{W^J^^^' = 
0.6) = 0.24], ^{W^fJ>^ = 0.3) = 0.78, and j{W^f:'' = 
0.6) = 1.02. The results of equation (2) for W^fJ^^ = 0.3 
A and 0.6 A arc shown in Figure 28. For comparison, SS92 
show in their figure 4 the number of sight-lines for these 
same values of W^^J^^; this is equivalent to redshift path 
density. 

Using equation (1), the effective redshift path for the 
RT survey is found to be AZnTiWrmn = 0.3) = 104.6 
with mean redshift < z >rt= 0.83 and 87 systems, and 
AZiiT{Wmin = 0.6) = 103.7 with the same mean redshift 
< z >RT= 0.83 and 44 systems. To confirm these re- 
sults we have verified that the detected number of Mg II 



^fMgII{W, 



A2796 
min 



Jo 



mRT{z,W^il'^)dz 



f 

Jo 



9RT{z,W^Z'')n{z,W^t'^)dz 



A2796\ 



(4) 



where AfMgii{W^T' = 0.3) = 87 and A4f,7/(W^™ = 
0.6) = 44. For comparison, we give the number of sys- 
tems, mean redshifts, and redshift paths of the RT and 
SS92 samples in Table 6. The fact that SS92 were able 
to study w 28% more systems than in the RT survey over 
a redshift path length that was only w 9% larger is un- 
derstandable since the SS92 survey had < z >ss92= 1.12, 
while the RT survey has < z >rt= 0.83. After correcting 
for the non-equivalent redshift paths for different W^^J^^ 
thresholds, we find that (per unit redshift) 55% of the 
Mg II systems in SS92 have Wo^2796> q q ^Ydle this frac- 
tion is « 50% for the RT survey. This is consistent with 
expectations since the SS92 and RT samples are derived 
from the same parent population of absorbers (§3.1). 

These determinations should be useful in attempts to 
compare the effectiveness of our survey strategy, and our 
results, to conventional surveys. At the present time, the 
observed incidence of DLA in the various surveys appears 
to be acceptably consistent given the small number statis- 
tics (Jannuzi et al. 1998). However, to make a reliable pre- 
diction of the number of systems which will be detected in 
a survey one should explicitly make use of equation (4) in- 
stead of simply multiplying the total redshift path, AZ, at 
a mean redshift < z > by the expected number of systems 
per unit redshift, n(< z >). 

4. THE LOW REDSHIFT DLA STATISTICS 

In this section, we first discuss the steps taken to cor- 
rect for any Malmquist bias in the RT survey (§4.1). We 
then derive results on the incidence of low-rcdshift DLA 
systems, noLA{z) (§4.2), and their corresponding cosmo- 
logical neutral gas mass density, i^oLAiz) (§4.3). We use 
the riDLAiz) and ^dla{z) notation deliberately to refer 
only to classical DLA systems having neutral hydrogen col- 
mnn densities Nhi > 2 x 10^" atoms cm^^. In particular, 
^DLA explicitly excludes hydrogen in molecular and ion- 
ized forms. It also excludes neutral hydrogen columns that 
may have been missed due to selection effects, such as ob- 
scuration of background QSOs by dust. Since Mg II lines 
are present in these systems, the gas is clearly enriched 
with metals. However, there is obviously insufficient dust 
to cause them to be overlooked. The presence of heavy 
elements in the neutral components of DLA systems does 
not make a significant contribution to ^Idla{z), however 
the presence of helium does. Therefore, helium from big 
bang nucleosynthesis is explicitly included in the calcula- 
tion of flDLA{z)- We assume a cosmology with A = and 
use Hq = 65 km s~^ Mpc~^ (has = Hq/65) for all cal- 
culations. Where applicable, we report results assuming 
qo = and qo = 0.5. 

4.1. Correction for Malmquist Bias 

A Malmquist bias can be introduced into a statistical 
study when a truncation or censoring effect infiuences sam- 
ple selection. This is a potential problem in absorption-line 
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Fig. 27.— The observed redshift distribution of Mg II systems in the RT sample, mRT{z, W^^^^), for W^y^^ = 0.3 A and 0.6 A. 
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Fig. 28.— The observed effective redshift path density of the RT survey, gRT{z,W^y^^), for W^T^^ = 0.3 A and 0.6 A. 



Table 6 

Number of Systems, Mean Redshift, and Redshift Path for the RT 

AND SS92 Samples 



Wmin (A) 




RT 






SS92 




NmqII 


< z> 


AZ 


NmqII 


<z> 


AZ 


0.3 


87 


0.83 


104.6 


111 


1.12 


114.2 


O.G 


44 


().<S3 


103.7 




1.17 


129 
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surveys since sample selection is usually based on a min- 
imum rest equivalent width or, for cases in which Voigt 
damping profiles are fitted to the data, the H I column 
density. This effect is significant if the column density is a 
steep-enough function near the threshold column density. 
For example, a Malmquist bias would be introduced if we 
identified and used only those systems with fitted Lya 
profiles that indicated Nhi > 2 x 10^° atoms cm~^ be- 
cause the measurements have associated errors and in re- 
ality some of the systems are likely to have Nhi < 2 x 10^° 
atoms cm~^. At the same time, there might be many more 
systems formally measured to lie below the H I column 
density threshold if errors were not taken into account, but 
in reality some of these are likely to lie above the threshold. 

Thus, our initial intent was to form a DLA sample which 
included those Mg II systems with measured H I column 
densities as small as 2cr less than 2 x 10^'' atoms cm~^. 
However, we have not identified and measured any systems 
which are slightly below the threshold. This is likely to be 
an artifact of small number statistics'\ Consequently, in 
order to approximately correct our sample for the effect of 
the Malmquist bias, we have simply removed half of the 
systems with Nhi ~ 2 x 10^° atoms cm~^, reasoning that 
about half of them are likely to have actual H I column 
densities slightly below the threshold as suggested by the 
measurement errors. The details are given below. 

Specifically, our study has yielded 87 Mg II systems for 
which we have information on the nature of the Lya line. 
Twelve of these systems were found to have DLA with col- 
umn densities Nhi > 2 x 10^° cm~^ and two more were 
found serendipitously (§2.5). However, as can be seen from 
Table 3, four systems have measurement errors that sug- 
gest that their true H I column densities may lie below 
the threshold value. These are the Nhi = 2.0 ±0.2 x lO^o 
atoms cm-2 z = 0.518 system toward B2 0827-1-24, the 
Nhi = 2.1 ± 0.5 x 10^° atoms cm'^ z = 1.391 sys- 
tem toward Q0957-h561A, the Nhi = 2.0 ± 1.0 x 10^° 



detections arise in the T4^o^2796> 0.6 A subsample (§3.1), 
so we believe that the effects of truncation at some Wq^''^^^ 
threshold are understood. Moreover, when we considered 
effects due to W^'^''^^ selection in RTB95, we found that 
when measurement errors were well-determined, the re- 
siilting statistical results remained unchanged for samples 
with different threshold 14^(^^2796 Therefore, we will ignore 
the effect of any Malmquist bias in the Mg II samples. 



4.2. UDLAiz) 

In §3.1 we showed that the low-redshift DLA absorbers 
are almost exclusively drawn from the population of ab- 
sorbers with Wo'^2^^^> 0.6 A. Therefore, we will use the 
y;j/A2796> 0.6 A Mg II Sample to derive the statistical prop- 
erties of DLA systems at low redshift. We start with the 
redshift evolution of the number density distribution of 
Mg II systems with Wo^2796> q.S A, nMgii{z,W^;i^^ = 
0.6). SS92 find that 



riMgiiiz, W^?™^ = 0.6) = (0.24 ± 0.10) (1 + z) 



1.02±0.53 



atoms cm z — 0.633 system toward Q1209+107. and 
the Nhi = 2.3 ± 0.4 x lO^'" atoms cm'^ z = 0.656 system 
toward 3C 336 (1622-1-239). As noted above, we found 
no systems which lie slightly below the H I column den- 
sity threshold. In order to correct for the Malmquist bias 
we have therefore excluded half of these four systems by 
giving them weight 0.5 when forming the DLA statisti- 
cal sample for the purpose of determining noLAiz) and 
^dla{z)- The correction will have the net effect of low- 
ering our UDLAiz) determination by 20% (§4.2), but in 
fact it has a minor effect on the determination of ^dla{z) 
since the contributions to flDLA{z) come from the highest 
column density systems. After correcting for Malmquist 
bias, we effectively have 10 DLA systems among 87 Mg II 
systems with Wq'^'^^^> 0.3 A and 9 DLA systems among 
44 Mg II systems with Wo^2796> q g ^ 

It should be noted that, in principle, a Malmquist bias 
must also affect the Mg II samples themselves since they 
are defined using threshold values for Wq"^"^^^. Since 
our particular Mg II sample was compiled from various 
sources, it is difficult to assess the degree to which the 
Malmquist bias affects the sample selection for Wq'^'^^^> 
0.3 A. However, with the exception of one system, our DLA 

^Based on the column density distribution we derive in §5, the expected 
atoms cm~^ is ~ 2. The probability of finding none is 13.5%. 



(5) 

nDLA{z) can then be determined by noting the fraction 
of DLA absorbers among the Mg II absorber population. 
Given the small number of DLA systems, we can reason- 
ably split the redshift range, 0.11 < z < 1.65, into at 
most two intervals. Recall that z = 0.11 is the low red- 
shift limit of optical spectroscopic surveys for Mg II and 
z = 1.65 is approximately the upper redshift limit of UV 
surveys for Lya. The mean redshift of the 44 Mg II sys- 
tems with Wo^2796> 0.6 A is < 2 >= 0.83. We observed 
21 Mg II systems in the interval 0.11 < z < 0.83 with 
< z >= 0.49. Employing the half-weighting method to 
correct for Malmquist bias (§4.1), three DLA systems in 
this redshift interval are assigned weight 1.0 and three oth- 
ers have weight 0.5. We observed 23 Mg II systems in the 
interval 0.83 < z < 1.65 with < z >= 1.15. Four of these 
are DLA systems with an assigned weight of 1.0, and one 
other has weight 0.5. Thus, using equation (5) we find. 



riDLAiz = 0.49) = ( 



,3 



-2.92 
-1.66 



-hO.5- 



h2.92 
-1.66 ^ 



21 21 ' 

(0.24 ± 0.10)(1 -I- o.49)^-°2=^°-^3 
- 08+° °^ 



(6) 



and 



riDLAiz = 1.15) = ( 



.-1-3.16 
*-1.91 



1 



0.5- 



+2.30 
-0.83 



23 23 
(0.24 ± 0.10)(1 -I- 1.15)^-°^=^°-^^ 

= 0.10t°:jo (7) 

where errors derived from small-number Poisson statistics 
are reported (Gehrels 1986). These results are shown in 

Figure 29. 

Figure 29 shows the overall results on riDLAiz) in the 
redshift interval < z < 4.7. The logarithm of nDLA{z) 
is plotted against log(l -|- z). The corresponding linear 
redshift scale is also shown. The points at z > 1.65 
are taken directly from WLFC95 and Storrie-Lombardi 
& Wolfe (2000, in preparation, kindly communicated to 

number of systems with column density 1.6 x 10^" < Nui < 2 x 10^" 
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us by L. Storric-Lombardi) . The highest-redshift data in- 
clude recent confirmed DLA detections. The z = point is 
derived from the observed H I distribution in local spiral 
galaxies (RTB95) and represents the probability of inter- 
cepting H I columns with Nhi > 2 x 10^" atoms cm~^. 

Due to small number statistics, the error bars shown in 
Figure 29 are significant, exceeding approximately 60% in 
the low-redshift bins and approximately 20% in the mod- 
erate to high redshift bins (excluding the highest redshift 
bins). Nevertheless, it is appropriate to address the issue 
of whether the new low-redshift result on nDLA{z) can be 
directly compared to the higher redshift results as done 
in Figure 29. This is because, even after setting aside 
the issue of possible biases (e.g. biases due to gravita- 
tional lensing or dust obscuration) , the completeness of the 
higher redshift surveys should be considered. These sur- 
veys include the ones described in WTSC86, Lanzetta et 
al. (1991), LWT95, and WLFC95. WLFC95 explain that 
their results are derived from the statistical analysis of 62 
DLA systems, but that "the complete statistical sample 
is likely to exceed 62 DLA systems, with additional sys- 
tems coming from the 102 untested candidates with (Lya 
absorption rest equivalent widths) 5 < W < 10 A." They 
further explain that "the 102 untested candidates arc likely 
to be confirmed (to be classical DLA absorption systems) 
at a rate lower than one in four." Thus, the problem with 
directly comparing the low-redshift results on nDLA{z) to 
the high- redshift residts (Figure 29), is the absence of a 
correction for the high-redshift incompleteness. One way 
to account for high-redshift incompleteness is to perform a 
statistical analysis which uses the considerably more com- 
plete subset of DLA systems which have Nhi > 4 x 10^° 
atoms cni^^ (e.g. as done in Lanzetta ct al. 1991). How- 
ever, since our low-redshift survey does not suffer from this 
problem, the method we adopt here is to correct for the 
high-redshift incompleteness. The details of this correction 
are discussed below and the corrected noLAiz) results are 
shown in Figure 30. 

First, we note that eight of the 102 untested candidates 
referred to by WLFC95 come from the low-redshift lUE 
sample of LWT95. In addition, four of the 62 DLA systems 
in the WLFC95 statistical sample have redshift z < 1.65. 
This suggests that there are 58 confirmed high-redshift 
DLA systems in the WLFC95 statistical sample plus an- 
other 94 "untested" candidates as of 1995. Our current 
count of untested high-redshift candidates shows that 89 
remain. To assess how the untested candidates may af- 
fect the high-redshift nDLA{z) results, we note the fol- 
lowing: (1) in the WLFC95 1.5 < z < 2.0 bin there are 
18 confirmed DLA systems and 33 untested candidates 
(in reality this bin can only extend down to z « 1.65), 
(2) in the WLFC95 2.0 < z < 2.5 bin there are 25 con- 
firmed DLA systems and 39 untested candidates, (3) in 
the WLFC95 2.5 < z < 3.0 bin there are 11 confirmed 
DLA systems and 12 untested candidates, and (4) in the 
WLFC95 3.0 < z < 3.5 bin there is one confirmed DLA 
system and 5 untested candidates. If the effective success 
rate for untested candidates (i.e. minus the correction for 
"confirmed" DLA systems that turn out to lie slightly be- 
low the Nhi — 2x 10^° atoms cm^^ threshold) is taken to 
be 20%, this indicates that the high-redshift nuLAiz) re- 
sults in WLFC95 should be increased significantly relative 



to the reported riDLAiz) errors (sec Figure 29). Evidence 
for such a high effective success rate is readily apparent in 
figure 2 of Lanzetta et al. (1991). Thus, we find correc- 
tion factors of 1.37 ± 0.14, 1.31 ± 0.11, and 1.22 ± 0.14 for 
bins (1) — (3), respectively. The errors on the correction 
factors suggest that differences between the correction fac- 
tors for different bins are not significant. Thus, an overall 
correction factor of 1.31 ± 0.07 appears to apply to the 
high redshift data of WLFC95. The new corrected re- 
sults on n£)LA{z) in the redshift interval < z < 3 are 
shown in Figure 30 and we adopt them for the remainder 
of our analysis. Results for z > 3 are not considered since 
the high-redshift study of Storrie-Lombardi & Wolfe is in 
preparation and the correction factor, if required, cannot 
be assessed. This new result is still qualitatively consistent 
with a linear decline in the incidence of DLA systems with 
decreasing redshift. 

While the statistics are poor, trends can nevertheless 
be extracted from the data. The decline in noLAiz) can 
be expressed using the parameterization that is normally 
adopted to describe the statistics of QSO absorption-line 
systems, 

nDLA{z)=no{l + zy. (8) 

Therefore, several forms of this power-law parameteriza- 
tion are shown on Figure 30. First, a single power-law is 
overplotted to include all of the data, including the z = 
point (i.e. the open circle, RTB95). The power- law has 
index 7 = 2.5, but a single power law may not fit all of the 
data satisfactorily. Second is a single 7 = 1.5 power law 
overplotted on all of the data except the z = point. This 
very adequately describes the incidence of DLA systems in 
QSO absorption-line surveys, but the extrapolation over- 
predicts the inferred incidence of DLA systems at z = by 
a factor of » 3. Third is a broken power-law overplotted on 
all the data, including the z = point. For convenience, 
we have assumed that the break occurs at redshift z = 1.5, 
consistent with the determination of the position of the 
break from the HST QSO Absorption-Line Key Project 
data which applies to the lower column density Lya forest 
lines (Weymann et al. 1998). However, at z < 1.5 the Key 
Project determined 7 « 0.16 for the Lya forest but we find 
that a much steeper decline, with 7 w 3.0, applies for the 
DLA systems if the z = point is included. At z > 1.5 
Bechtold (1994) found 7 w 1.85 for the Lya forest but we 
find 7 « 0.4 for the DLA systems. Thus, the data indicate 
that the broken power-law indices for the DLA systems are 
significantly different from the indices determined for the 
Lya forest. For the DLA systems, the incidence is nearly 
fiat at moderate-to-high redshift, but then must decline 
rapidly at low redshift to fit the z = point, while the 
Lya forest shows a more rapid decline at high redshift but 
then flattens at low redshift. However, the Key Project 
data do show evidence for a steepening in the incidence 
of Lya forest lines with increasing rest equivalent width. 
The incidence of DLA systems at z < 1.5 is qualitatively 
consistent with this trend, but it is not consistent with 
this trend at z > 1.5. Moreover, any break in the DLA 
power-law fit may actually occur at some redshift other 
than 1.5, but this cannot be determined from the present 
data. 

What do these results imply? Recall that for A = 
cosmologies, the comoving number density of absorbers is 
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Fig. 29. — A plot of the logarithm of noLAiz), the number of DLA systems per unit redshift along a line of sight, as a function of the 
logarithm of {1 + z). The corresponding linear redshift scale is also shown. The new low-redshift results are shown as solid squares. The solid 
circles are from WLFC95 and the solid stars are from Storrie-Lombardi &; Wolfe (in preparation) . These data include recent confirmed DLA 
detections. The 2 = point is derived from the observed H I distribution in local spiral galaxies {RTB95) and represents the probability of 
intercepting H I columns with Nhi > 2 x 10^" atoms cm~^. Vertical bars correspond to Icr errors and horizontal bars indicate bin sizes. 




log(l + z) 

Fig. 30. — Same as Figure 29 but with a Malmquist bias correction factor applied to the high-icdsliift data of WLFC95 (sec text). Results 
for z > 3 arc not included since the high-rcdsliift study of Storrie-Lombardi & Wolfe is in preparation, and the correction factors, if required, 
cannot be assessed. Three different forms of the power-law fit are shown (see equation (8)). The solid line has 7 = 2.5 and is forced to pass 
through the z = data point, the dotted line has 7 = 1.5 and does not include the z = data point. Extrapolation of this power-law to 
z = (the open triangle) results in a value that is ~ 3 times larger than the observed incidence at z = 0. The dashed line has 7 = 0.4 for 
z > 1.5 and 7 = 3.0 for z < 1.5. 



32 



consistent with no intrinsic evolution if 7 = 0.5 for qo = 0.5 
or if 7 = 1.0 for qQ = 0. Thus, if we neglect the z = point, 
the present data on nDLA^z) with 7 = 1.5 shows some ev- 
idence for real intrinsic evolution, but the statistics are 
not good enough to definitively rule out no evolution in 
a go = universe. However, if a no-evolution trend is 
valid down to z = 0, it implies that the H I 21 cm emis- 
sion surveys are missing « 70% of gas-rich systems with 
Nhi > 2 X 10^*^ atoms cm~^ locally. We consider this pos- 
sibility unlikely since HI 21 cm emission surveys have not 
revealed a significant population of gas-rich objects that 
are not already included in optical surveys (Briggs 1997; 
Zwaan et al. 1997). These H I surveys are sensitive down 
to Nhi > 10^^ atoms cm~^ and are capable of detecting 
masses Mhi > 2 x Kfh^^ M,©. Alternatively, the need 
for a large value for the power-law index, e.g., 7 ~ 2.5 or 
the broken power law with 7 = 0.4 and 7 = 3 (2: > 1.5 
and z < 1.5, respectively), arises solely because of the in- 
clusion of the z = data point. Both this steep single 
power-law or the broken power-law imply that the DLA 
absorbers are undergoing significant evolution. At z « 2 
the incidence of DLA systems is a factor « 17 times larger 
than inferred for z = 0. This is « 6 times larger than 
what is predicted in a go = no-evolution universe (« 10 
times larger if go = 0.5). Since nDLA{z) is proportional to 
the number of absorbers per comoving volume times their 
H I cross-section, the decline in incidence suggests that at 
z < 2, during the last 67% (go = 0) of the age of the Uni- 
verse (81% for go = 0.5), there was either a decrease in the 
effective H I cross-sections of DLA absorbers, a decrease 
in the number of DLA absorbers per comoving volume, or 
a combination of both effects. A decrease in cross-section 
may correspond to a collapse phase for DLA absorbers, 
while a decrease in the number of DLA absorbers may be 
the result of mergers among them. 

4.3. ftDLAiz) 

The cosmological mass density of neutral gas in DLA 
systems, iloLAiz), is determined using nDLA{z) and the 
mean column density of DLA systems as a function of red- 
shift, {Nhi{z)), as follows 



^dla{z) 



IJ.87TGmH nDLA(z){NHi{z)) 
3cHo (1 + z) 



for go = and 



DLA 



_ fj.8TrGm,H nDLA{z){NHi{z)) 
~ 3cHo (1 + z)i/2 



(9) 



(10) 



for go = 0.5. Here fj, = 1.3 corrects for a neutral gas 
composition of 75% H and 25% He by mass, G is the grav- 
itational constant, and mn is the mass of the hydrogen 
atom. A Hubble constant of Hq = 65 km s^^ Mpc^^ is 
used for all calcidations. To determine ClDLAiz) at low 
redshift we use the DLA systems and associated column 
densities given in Table 3. As before, we split the redshift 
interval 0.11 < z < 1.65 into two redshift bins. Three 
DLA systems with weight 1.0 and three DLA systems with 
weight 0.5 are in the lower redshift bin; four DLA systems 
with weight 1.0 and one with weight 0.5 are in the higher 
redshift bin. For go = 0, we find 



and 

^dla{z = 1.15) = {lA5t{f^) X 10-3. 
For go = 0.5, we find 

^dla{z = 0.49) = (2.23tiJ^) x lO'^ 



and 



ftoLAiz = 1.15) = {2.13tl,\l) X 10" 



\dla{z = 0.49) = (l.82tl;^8) x 10"=^ 



The logarithm of Qdla{z) for go = and go = 0.5 is 
plotted in Figure 31. The high-redshift results of WLFC95 
and Storrie-Lombardi & Wolfe (2000, in preparation) are 
shown as solid circles and stars, respectively. In this case, 
the completeness correction (see §4.2) to the high-redshift 
results on ^dla{z) is small (« 7%) relative to the reported 
errors and has not been applied to the data. See §6 for a 
discussion of what are likely to be more significant errors 
in the determination of ^l]JLAiz). The new low-redshift 
points are plotted as solid squares and are seen to be of 
comparable value to the points at higher redshift. The er- 
rors on the low-redshift points are large owing to the small 
number of low-redshift DLA systems studied. The open 
circle at z = is the local neutral gas mass density as mea- 
sured by Rao & Briggs (1993). The open square at z = 
is the local luminous mass density in stars as determined 
by Fukugita, Hogan, & Peebles (1998). 

The results on il]jLA{z) are similar to the nDLA{z) re- 
sults in that they show strong evidence for evolution if the 
z = point is included. If go = 0, the comoving cosmo- 
logical mass density of high-redshift DLA absorbers is « 4 
times larger than what is inferred from the local 21 cm ob- 
servations (« 6.5 times larger if go = 0.5). Alternatively, 
if the z = point is excluded, there is no convincing evi- 
dence for a decrease in flDLA{z) with decreasing redshift 
even at redshifts as low as z w 0.5. In fact, the actual 
results reveal a tendency for ^lDLA{z) to rise slightly with 
decreasing redshift in the interval z ~ 2—0.5. As expected, 
this tendency is somewhat more pronounced in a go = 
universe than in a go = 0.5 universe. However, the error 
bars on the Qdla{z) low-redshift points are fractionally 
much larger than the noLAiz) error bars and therefore 
more problematic. This is because the determination of 
^dla{z) is very sensitive to the small number of systems 
which have the highest column densities. Thus, the cur- 
rent results do not provide very strong constraints on the 
evolution of the amount of neutral gas in the Universe in 
the redshift interval 0.1 < z < 1.6. 

Nevertheless, assuming that the z = point is correct 
(i.e. assuming we have not missed a significant fraction of 
the z = neutral gas in local 21 cm emission surveys), the 
data would be consistent with the possibility of a rapid 
decline in both nDLA{z) and i^DLAiz) from z « 0.5 to 
z = 0. If go = 0, this would imply that the DLA absorber 
population lost w 75% of its neutral gas mass over the 
most recent w 33% of the age of the Universe (if go = 0.5, 
« 85% of the mass would be lost during the most recent « 
46% of the age of the Universe). We also note that, except 
for the z = 4.1 data point, all the z > data points in 
Figure 31 lie within la of constant logflDLAiz) = —3.0 
for go = and logiloLAiz) = —2.8 for go = 0.5. The 
z = 4.1 data point is discrepant from a constant value at 
about the 2a level. 
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Fig. 31. — A plot of the logarithm of ^dla{z) with Ho = 65 km Mpc^^ for qo = (top panel) and go = 0.5 (bottom panel). The new 
low-redshift results are shown as solid squares. The solid circles are from WLFC95 and the solid stars are from Storrie-Lombardi Sz Wolfe 
(in preparation). Vertical bars correspond to la errors and horizontal bars indicate bin sizes. The open circle at z = is the local neutral 
gas mass density as measured by Rao & Briggs (1993) and the open square at z = is the local luminous mass density in stars (Pukugita, 
Hogan, & Peebles 1998). 
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What do these results imply? In the context of the con- 
version of neutral gas into stars and global galaxy forma- 
tion as might be tracked by ilDLA{z), the current data 
provide no evidence which would support the view that the 
bulk of galaxy formation happened, for example, at red- 
shifts z > 1. This is in contradiction with earlier conclu- 
sions reached by LWT95 on low-redshift DLA; their re- 
sults showed a steady decline in ^dla{z) from 2; « 3.5 to 
z w 0.008. The constant neutral gas mass density from 
z ~ 4 to z ?a 0.5 that our new results imply is also at 
odds with the evolution of the star formation density in 
the Universe as derived from the UV luminosity density 
of galaxies. Steidel et al. (1999) show that after cor- 
rection for effects of dust obscuration, the star formation 
density remains roughly constant from z « 4 to z « 1. 
The data for z < 1 from Lilly et al. (1996) indicate that 
the star formation density falls by a factor of w 5 from 
z » 1 to z « 0.2. The star formation density and the 
gas mass density cannot both be constant over the same 
range of redshifts (1 < z < 4) if they track the same ob- 
jects. This implies that DLA galaxies and star- forming 
galaxies trace different populations (see also Pettini et al. 
1999, Steidel et al. 1999 and references therein). However, 
star-forming galaxies must have neutral gas from which to 
make stars. There are possibly three reasons why searches 
for DLA systems are not finding star-forming galaxies: (1) 
the intensive star-forming phase of UV-luminous galaxies 
is short-lived, (2) background quasars are obscured by dust 
in the UV-luminous galaxies, and/or (3) the UV-luminous 
galaxies are likely to be the brightest galaxies at the epochs 
they are detected (whereas, random sight lines that pick 
up the DLA population are more likely to intersect H I re- 
gions in the more common type of object, i.e., L* or fainter 
galaxies. Additional evidence from metallicity studies and 
direct imaging also point to DLA galaxies as being a slowly 
evolving population. We will return to the discussion on 
the nature of DLA galaxies in §7. 

5. THE COLUMN DENSITY DISTRIBUTION FUNCTION OF 
DLA SYSTEMS 

The column density distribution function, f{NHi,z) = 
d'^Af{z) I dNnidX , specifies the number of absorbers per 
unit H I column density per unit absorption distance, 
where the absorption distance is X{z) = |[(1 + z)^ — 1] 
if 90=0 and X{z) = |[(1 zf-^ - 1] if go=0.5. Thus, 
f{NHi , z)dNHidX is the number of absorbers with H I 
column density between Nhj and Nhi + dNni in absorp- 
tion distance interval dX. 

There are several reasons why it is useful to study this 
function. For example, its shape can reveal information 
on the geometry of the absorbing regions (e.g. Milgrom 
1988). More ideally, evolutionary models which also take 
into account the global physical state of the gas (the H+, 
H'', and H2 fractions) and the conversion of gas into stars 
would predict fiNni-z). This prediction could then be 
empirically tested. Unfortunately, the statistics are still 
too poor to realize this goal in detail. However, we can 
form a sample which will permit us to derive the most ac- 
curate results currently possible on the H I column density 
distribution function of low-redshift DLA systems (§5.1). 
We accomplish this by adopting a rationale that allows 
us to form the largest possible unbiased low-redshift sam- 
ple. After correcting the sample for Malmquist bias, the 



normalized ciimulative column density distribution can be 
derived. We then use the KS test to compare the low- 
redshift normalized cumulative column density distribu- 
tion function to ones derived from high-redshift DLA sys- 
tems (§5.2) and local spirals (§5.3). Most notably, we show 
that the comparisons reveal clear evidence for differences 
in the shapes of these three column density distribution 
functions, which we interpret as being indicative of a ma- 
jor evolutionary trend. These results are then converted 
to absolute DLA column density distribution functions to 
show the redshift evolution (§5.4). 

5.1. The Low-Redshift DLA Column Density Distribution 

5.1.1. Rationale for Sample Selection 

The selection criteria for constructing a DLA column 
density distribution are different than what were used to 
define the RT, and hence, the RT-DLA samples (§3). We 
now explain how this sample was constructed. To derive 
the most accurate results possible on the low-redshift DLA 
column density distribution function we should use all 
available unbiased information on all known low-redshift 
DLA systems. Table 7 lists all of the confirmed low- 
redshift (z < 1.65) classical DLA {Nhi > 2 x 10^° atoms 
cm~^) systems from cosmologically intervening neutral gas 
that are currently known to us, without regard to possi- 
ble bias. By confirmed systems we mean ones for which 
good UV data have been obtained on the Lya absorption- 
line profile or for which 21 cm absorption lines have been 
detected. For both of these cases, a reliable estimate of 
the H I column density can be made. This is preferably 
done by fitting a Voigt profile to the damping wings of the 
Lya absorption line (§2), although use of the radio data 
and the assumption that the gas has 21 cm spin tempera- 
ture of Ts w 700 K can also be used to produce a reliable 
estimate (Lane et al. 2000, in preparation). Thus, the 
systems included in Table 7 form a biased sample of low- 
redshift classical DLA systems. Excluded from Table 7 are 
any non-classical "damped Lya systems" with H I column 
densities more than 2cr less than 2 x 10^*^ atoms cm~^ (al- 
though damping wings on the Lya absorption-line profile 
may still be easily visible in such systems), and all poten- 
tial DLA systems which have been proposed on the basis 
of certain metal-line ratios, x-ray absorption, or still even 
less-direct evidence. Also excluded are 21 cm absorption 
systems that are not cosmologically intervening along the 
sight-line to the quasar, but instead more likely arise in 
the quasar's host galaxy (e.g. Carilli et al. 1998). 

To derive the correct shape for the column density dis- 
tribution function, all systems identified on the basis of 
a bias in H I column density must be excluded. In Ta- 
ble 7 we note which of the systems were selected on the 
basis of such a bias. The remaining systems in Table 7 
form our "unbiased sample." As before (§2.4), we ex- 
clude the four previously known 21 cm absorbers. We 
also exclude cosmologically intervening 21 cm absorbers 
that have been identified in red quasars (e.g. Carilli et al. 
1998 and references therein) since their selection is based 
on the bias of observing red (i.e. dust reddened) quasars. 
Three systems fit into this category, and in some cases 
these also involve gravitational lensing which might pro- 
duce an additional bias. On the other hand, over the last 
several years four low-redshift DLA systems were identi- 
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Table 7 

Confirmed Low-Redshift Classical DLA Systems 



QSO 


V 




zdla 


7Vh//1020 
atoms cm~^ 


Notes 


Reference 


m /I Q m 


K .1 


0.14 


i.Dio 


z.U 




1 


noi Q 1 QKi 
Uzio+Jt) 1 


on n 
/U.U 


u.y4 


U.D8 / 


1 '} 
16 


Biased (red quasar) 


o 
z 


n o o r 1 ^ a A 
Uz35+iD4 


15.5 


0.94 


0.5z4 


50 


rJiased (zl cmj 


o 

o 


0248+430 


17.6 


1.31 


0.394 


36 


1 


0302-223 


16.4 


1.41 


1.010 


2.3 




1 


0454+039 


16.5 


1.34 


0.859 


4.7 




1 


0500+019 


21.2 


>0.58 


0.585'^ 


43'" 


Biased (red quasar) 


4 


0738+313 


16.1 


0.63 


0.091 


15 




1,5 








0.221 


7.9 




1,5 


0809+483 


17.8 


0.87 


0.437 


6.3 


Biased (21 cm) 


6 


0827+243 


17.3 


0.94 


0.518 


2.0 


1 


0933+732 


17.3 


2.52 


1.478 


42 




1 


0935+417 


16.3 


1.98 


1.372 


3.3 




7 


0952+179 


17.2 


1.47 


0.239 


21 




1 


0957+561 


17.0 


1.41 


1.391 


2.1 




1 


1122-168^= 


16.5 


2.40 


0.681 


2.5 




8 


1127-145 


16.9 


1.19 


0.313 


51 




1 


1209+107 


17.8 


2.19 


0.633 


2.0 




1 


1229-021 


16.8 


1.04 


0.395 


5.6 


Biased (21 cm) 


6 


1328+307 


17.3 


0.85 


0.692 


15 


Biased (21 cm) 


6 


1354+258 


18.0 


2.00 


1.418 


32 


1 


1622+239 


17.5 


0.93 


0.656 


2.3 




1 


1830-211 


22.0 


0.89 


0.192 


2.5 


Biased (red quasar) 


9 



* This system could be the result of absorption by the quasar host galaxy. 

^ The column density was determined assuming a spin temperature of Tg = 700K (Lane et al. 2000, in preparation). 

de la Varga & Reimcrs (1997) also detect a DLA candidate towards Q0515— 442 at z = 1.15 using lUE spectra. 
In the past, lUE spectra have been found to be unreliable sources for DLA candidates and we have therefore not 
included this system in Table 7. 

References. - (1) This work, (2) CariUi, Rupcn, & Yanny 1993, (3) Cohen ct al. 1999, (4) CariUi ct al. 1998, (5) 
Rao & Turnshek 1998, (6) Boisse et al. 1998, (7) Jannuzi et al. 1998, (8) de la Varga & Reimers 1997, (9) Wiklind 
& Combes 1998 
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fied in the Lya forest sercndipitously, without regard to 
column density bias. The DLA systems discovered in this 
way were not identified using a technique which would bias 
the H I column densities to higher values. Therefore, these 
serendipitously-discovered systems are included in our un- 
biased sample. Specifically these arc the z — f .372 DLA 
system in Q0935+4f 7 which was found in the EST QSO 
Absorption Line Key Project survey (Jannuzi et al. 1998), 
the z = 0.681 DLA system in Q1122-168 (de la Varga & 
Rcimcrs 1997), and the z = 0.0912 system in 01 363 and 
the z — 1.613 system in UM 366 that we found in our HST 
spectra. (Note that these four systems do not have Mg II 
information and, therefore, are not part of the RT and RT- 
DLA samples.) In summary, the present total low-redshift 
sample consists of 23 low-rcdshift DLA systems (Tabic 7) , 
seven of which are excluded due to selection bias. Thus, 
before correction for Malmquist bias, we have amassed an 
unbiased sample of 16 systems to study the low-redshift 
DLA column density distribution function. 

The reasons for making a correction for Malmquist bias 
were discussed in §4.1. To approximately correct our sam- 
ple for this effect we have given the five systems with 
Nhi ~ 2 X 10^" atoms cm~^ a weight of 0.5 in the col- 
umn density distribution analysis. 

5.1.2. Normalized Cumulative Column Density 
Distribution Function 

Shown in Figure 32 as a solid curve with open circles is 

the normalized low-redshift DLA cumulative column den- 
sity distribution function. This was derived using the 
Malmquist-bias-corrected sample discussed above which 
has mean redshift < z >« 0.78. Since the curve is normal- 
ized, we can compare the shape of this curve to the shapes 
of similarly derived column density distribution functions 
for high-redshift DLA absorbers and local spirals. This 
will be done in §5.2 and §5.3, respectively. 

5.2. Comparison with Results at High Redshift 

In order to compare our low-redshift result with results 

obtained at higher redshift, we have derived the high- 
redshift DLA normalized cumulative column density dis- 
tribution function using the data of WLFC95. Wc did 
this by forming an unbiased sample of high-redshift col- 
umn densities as follows. First, we included all 43 high- 
redshift (z > 1.65) systems with higher-resolution follow- 
up spectroscopy which have been used to confirm, with lit- 
tle doubt, that these systems are damped. In these cases 
H I column densities were derived by fitting Voigt damp- 
ing profiles to the Lya absorption lines. We then included 
all 15 candidate DLA systems that have Lya rest equiva- 
lent widths > 10 A. These systems do not have follow- 
up spectroscopy, but their large rest equivalent widths 
formally indicate that they have Nhi > 2 x 10^° atoms 
cm~^, although a few may not. Wc converted these rest 
equivalent widths to column densities using the relation 
Nhi = 2(Wo^i2i6/ioA)i/2 x 10^° atoms cm'^. These 58 
systems comprise the high-redshift {z > 1.65) statistical 
sample of WLFC95. All but one of these systems lies in 
the redshift interval 1.65 < z < 3.0. Finally, in keep- 
ing with the procedure adopted in §4.2, we assumed that 
18 of the remaining 89 "untested" candidate systems (« 
20%) are damped and, in reality, have Nhi > 2 x 10^° 



atoms cm^^. When adding these 18 systems to the unbi- 
ased high-redshift sample, we have included them evenly 
over the column density interval 20.3 < \og{NHi) < 20.4, 
consistent with likely measurement errors. Recall that 
WLFC95 conclude that they expect these untested can- 
didates to be confirmed at a rate lower than one in four. 
However, evidence for an effective success of at least 20% 
is readily apparent in figure 2 of Lanzetta et al. (1991). 
We chose a success rate of 20% since some of the systems 
included in the first two steps may actually have Nhi 
slightly smaller than 2 x 10^" atoms cm~^. This makes 
an approximate correction for the Malmquist bias in the 
high-redshift sample. The resulting normalized cumulative 
column density distribution function derived from this un- 
biased high-redshift sample is plotted on Figure 32 along 
with the low-redshift result. The high-redshift result is 
derived from an effective sample size of 76 DLA systems 
with mean redshift < 2; >« 2.26. 

A KS test comparing the low- and high-redshift results 
indicates that there is only a 2.8% probability that the 
two samples are derived from the same parent popula- 
tion. This comparison refers only to the shapes of the col- 
umn density distribution functions and not their normal- 
izations. The low-redshift sample clearly contains a much 
larger fraction of high column density systems. Thus, our 
analysis indicates that the interpretation of the current 
data are inconsistent with an earlier interpretation of the 
H I column density distribution of low-redshift DLA sys- 
tems (LWT95). Further discussion of this result is deferred 
to §5.5. 

5.3. Comparison with Results Derived from Local Spiral 

Galaxies 

Rao & Briggs (1993) used H I emission observations of 
27 local spirals and information on the spiral galaxy lumi- 
nosity function to construct a local H I column density 
distribution function. The more recent study of Zwaan et 
al. (1999) is consistent with the Rao & Briggs result. A 
comparison of the Rao & Briggs (1993) result to the one 
derived from high-redshift DLA statistics by Lanzetta et 
al. (1991) showed clear evidence for evolution in the num- 
ber of systems with redshift, but the error bars on the high- 
redshift points did not permit any meaningful conclusions 
to be drawn from a comparison of the shapes of the local 
and high-redshift H I column density distributions. Now, 
with improved statistics from high-redshift DLA systems 
(WLFC95) and the results on low-redshift DLA systems 
reported here, it is possible to perform a more meaningful 
comparison of the shapes of the column density distribu- 
tion functions versus redshift. The one caveat which might 
influence comparisons with the local results is the fact that 
the beam size used in H I 21 cm emission surveys might 
be too large to resolve out the local H I cloud structure. 
If this effect were severe enough, it would cause the high- 
est column densities to be under-represented in the local 
21 cm emission surveys. We believe the importance of 
this effect is minimal and discuss it in more detail in §6.5. 
Of course, this is not a problem in the pencil-beam QSO 
absorption- line surveys. The result on the local H I cumu- 
lative column density distribution derived by Rao & Briggs 
(1993) is included on Figure 32. A KS test comparing the 
low-redshift DLA result with the local H I result indicates 
that there is only a 0.01% chance that the two samples 
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Fig. 32. — The normahzed cumulative column density distributions for the low-redshift DLA sample (solid curve with open circles), the 
WLFC95 high-redshift DLA sample (dotted curve) , and local galaxies (solid curve) . The WLFC95 and RT-DLA samples have been corrected 
for Malmquist bias. Data for local galaxies are from Rao & Briggs (1993). 



are derived from the same parent population^. Again, the 
comparison refers only to the shapes of the column den- 
sity distribution functions and not their normalizations. 
The comparison shows that the local H I sample contains 
a much smaller fraction of very high column density gas 
than is found in the low-redshift DLA sample. In fact, the 
deficit of high column density gas in the local H I sample 
is even more pronounced than in the high-redshift DLA 
sample (§5.2). 

5.4. The DLA Absolute Column Density Distribution 

Function versus Redshift 

The absolute column density distribution functions at 
z = 0, low redshift (< z >« 0.78), and high redshift 
(< z >» 2.26) are shown in Figure 33 assuming qq = 
and go = 0.5. For z = this function is taken directly 
from Rao & Briggs (1993). The data points at low and 
high redshift are determined by first estimating the frac- 
tion of systems detected as a function of column density, 
y{NHi , z), for the low- and high-redshift samples, respec- 
tively. This is shown in Figure 34. Note the higher fraction 
of the highest column density systems in the low-redshift 
sample. Then, /{Nhj, z) is given by 



Thus, 



f{N. 



HI, 



= C{z. 



^ y{NHi,z) 

AN, 



HI 



(11) 



where the constant, C{z), is inversely proportional to the 
absorption distance interval of the survey (see §5). In turn, 

C{z) can be determined from noLA^z) at < z >= 0.78 and 
< z >= 2.26 by using the expression relating nDLA{z) and 

f{NHI,z), 



riDLAiz) = (1 + z){l + 2qoz) 



-V2^/(iV. 



hi,z)ANhi- 
(12) 



Ciz) = 



nDLA{z) 



(l + z)(l + 2go^)-i/2 Y^yi^Hi^z) 



(13) 



Note that ^ y(NHi, z) is just the total number of systems 
found in the survey. The resulting values of C{z) for the 
WLFC95 and RT samples are given in Table 8 for cos- 
mologies with go = and go = 0.5; /{Nhj, z) then follows 
from equation (11) and Figure 34. 

In Table 8 we also present the results of parametric fits 
to the data on the absolute column density distribution 
functions. These fits are over-plotted on the binned data 
(Figure 33). Given the statistical accuracy of the data, the 
results at both high and low redshifts are fitted by a single 
power law. It has been argued that such power-law fits 
are unphysical (e.g. WLFC95, Storrie-Lombardi, Irwin, & 
McMahon 1996a, Pei, Fall, & Hauser 1999) because, for 
example, they cause the mass in DLA systems to diverge 
when they take the form f{NHi,z) ~ N'^j with (3 < 2. 
While we agree that the power-law fits can be unphysi- 
cal, we have nevertheless opted to use them because the 
data at high and low redshifts do not offer support for 
any particular model. For example, the exponential-disk 
model fits to the column density distributions presented 
by WLFC95 are arbitrary in that the data provide no ev- 
idence for the breaks that are shown in their figure 7. In 
addition, at high redshift the data are in excellent agree- 
ment with a single power-law fit, while too few systems are 
known at low redshift to adopt a form other than a single 
power law. Moreover, as we discuss in §5.5, the shapes of 
the column density distribution functions at high and low 
redshifts appear to be inconsistent with a set of randomly 
distributed gas-rich disks. Therefore, we believe that it is 



®A KS test comparing the high-redshift DLA result with the local H I result indicates that there is only a 2.0% chance that these two 
samples are derived from the same parent population. 



38 



o 



qg = ■ This work 

« WLFC95 

^ ^ Rao & Briggs (1993) 




I , , , , 1 , , , , 1 , , , , 1 , , , , 1 

1 20 20.5 21 21.5 22 

log Njjj [cm^^] 

Fic. 33. — The logarithm of the absolute column density distribution function, f{NHi), for the low-redshift DLA sample (solid squares), 
the WLFC95 high-redshift DLA sample (solid circles), and local galaxies (solid line) for go = and qo = 0.5. Vertical bars correspond to Icr 
errors and horizontal bars indicate bin sizes. The dashed line is a power-law fit to the WLFC95 data with fS = 1.8 and the dotted line is a 
power-law fit to the RT data with /3 = 1.3. 
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Fig. 34. — The fraction of DLA systems detected as a function of column density. The bin size is 0.3 in logNuj. Note the higher fraction 
of the highest column density systems in the low-redshift sample. 
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Table 8 

Column Density Distribution Function Parameters 



Sample 


< z> 


90 = 


C{z) 

90 = 0.5 




WLFC95 


2.26 


(1.1 ±0.2) X 10" 


-3 (2.0 ± 0.4) X 10-3 


1.8 


RT 


0.78 


(3.9 ± 2.4) X 10" 


-3 (5.2 ± 3.2) X 10-3 


1.3 


Local galaxies 





1.2 (log Nhi < 21) 










2.9 (log Njjj > 21) 



more appropriate to parameterize the column density dis- 
tribution functions in a mo del- independent way. At both 
high and low redshifts, we adopt the view that there is a 
cutoff in H I column density at Nhi ~ 5 x 10^^ atoms 
cm-^. So far. this is the maximum H I column density 
found for a DLA system at any redshift. 

On the other hand, the absolute H I column density 
distribution at z = is clearly not wcU-rcprcscntcd by a 
single power law. There is clear evidence for a break in the 
column density distribution at Nhi ~ 10^^ atoms cm-^. 
We represent the data at 2; = using a broken power-law 
(Table 8), but again we assume that there is a cutoff in 
the H I column density distribution at Nhi » 5 x 10^^ 
atoms cm-^. 

To summarize, our results suggest that at Nhi ~ 2.2 x 
10^1 atoms cm-2, /{Nhi, < z >^ 2.2) w /{Nhi, < z >fv 
0.8); however /3 « 1.8 at < 2 >« 2.2, while P « 1.3 at 
< z >w 0.8. The value of f{NHi,z = 0) is correspond- 
ingly lower at all H I column densities, with /3 = 1.2 for 
Nhi < 10^^ atoms cm-^ and /3 = 2.9 for Nhi > 10^^ 
atoms cm-^. 

5.5. Discussion 

Despite our findings on the evolution of nDLA{z) and 
^dla{z), our results on the normalized DLA column den- 
sity distribution fimction, f{NHi,z), indicate that there 
is significant evidence for real evolution in /{Nhi, z) even 
after excluding the z = data points (Figures 32 and 
33). When we compared these trends using the normalized 
shape of the column density distribution function we found 
them to be highly significant (§5.2 and §5.3)^. For ex- 
ample, a comparison of the low-redshift and high-redshift 
normalized column density distribution functions indicates 
that there is only an « 2.8% chance that they are drawn 
from a population of absorbers with the same properties. 

This result may seem surprising since there was only 
marginal evidence for real evolution in nohAiz) and 
^dla{.z) per comoving volume (90 = universe) after the 
= data are excluded. However, a close examination of 
the data shows that the n^i^A^z) and ^£,i^a{z) results are 
consistent with the f(NHi,z) results. For a non-evolving 
9o = universe our results show that nDLAiz){l + z)~^ 
marginally decreases with decreasing redshift (i.e. in Fig- 
ure 30 a 7 = 1.5 power law fits the DLA data better than 
a 7 = 1.0 power law), while ^lDLA{z) marginally increases 
with decreasing redshift (Figure 31). These trends ulti- 

^Also, Storrie-Lombardi et al. (1996a) report some evidence for both 
the DLA incidence at the very highest redshifts (z ~ 4) that are reacha 



mately stem from the differences in /{Nhi, z) at high and 
low redshifts, and as discussed above, our results (Figure 
32 and 33) indicate that real evolution has occurred in 
fiNnitz). The results on nDLA{z) and 0£ila(z) do not 
explicitly contain any information as a function of column 
density. Thus, /{Nhi, z) provides new information which 
nDLA{z) and Q.dla{z) do not provide. 

What are the possible origins for the observed evolution- 
ary trend in the column density distribution Junction? We 
consider this question by discussing two conceptual mod- 
els for DLA absorbers. The first model is the Giant Hy- 
drogen Cloud (GHC) model proposed by Khersonsky & 
Turnshek (1996). In this model DLA galaxies are sim- 
ply presumed to contain GHCs that can be approximated 
by spheroids. The second model is the exponential-disk 
(ED) model which was recently discussed by WLFC95. In 
the ED model the H I column density falls off exponen- 
tially with increasing radius, in approximate agreement 
with 21 cm observations of H I in gas-rich spirals (e.g. 
van Gorkom 1993). Moreover, at increasing H I column 
densities, the column density distribution function must 
fall off as ~ ^Hi for randomly inclined disks (Milgrom 
1988; WLFC95). This is consistent with what is observed 
locally (Rao & Briggs 1993; Zwaan et al. 1999), where 
the /{Nhi) ~ N^^j behavior applies at Nhi > 10^^ atoms 
cm^^. In the ED model of WLFC95, it is supposed that 
the /{Nhi) ~ N^^^ behavior sets in at increasingly larger 
H I values with increasing redshift. We note that the 
GHC and ED models are not necessarily inconsistent with 
each other since it is possible to envisage roughly spherical 
GHCs arranged in a galactic disk. 

We will consider the case where the z = Q data must fit 
into any evolutionary trend, as this seems to be a reason- 
able constraint. 

What is the cause o/ the /lattening o/ the /{Nhi, z) dis- 
tribution /rom high to low redshift? First, recall that the 
observed flattening from < z 2.2 to < 2; >« 0.8 is just 
the opposite of the evolutionary trend reported by LWT95. 
We will consider the evolution in the context of both the 
GHC and ED models. For the GHC model, one would ex- 
pect /{Nhi, z) to flatten if accretion processes were taking 
place that caused the relative build-up of more massive H I 
clouds in the central regions of the absorbing structure. As 
an illustrative example, consider a case noted by Milgrom 
(1988). Namely, he points out that for random sight-lines 

a steepening in the column density distribution function and a decrease 
ble in surveys. 
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through identical spherical systems where the column den- 
sity at impact parameter I is Nhj(I) = A^o[l + (V^c)^]"*^: a 
power-law column density distribution with (3 = (c + l)/c 
would result for lines of sight that are outside the core. 
Here, A^o is the column density at ^ = 0, Tc is the core ra- 
dius, and c is the power-law index which specifies the run 
of the neutral gas density distribution with cloud radius, 
no{r) ~ [1 + {r/rc)^]~^'^~^^^^\ In this illustrative example 
we would expect c « 1.3 at < ^; >w 2.2 and c w 3.3 at 
< z >« 0.8. This corresponds to much more centrally 
condensed systems at lower redshift. Thus, in the con- 
text of this example, such processes might be associated 
with galaxy collapse or mergers of galaxies. In particu- 
lar, if this type of process corresponded to both a relative 
increase in the masses of GHCs and an increase in their 
particle density (i.e. somewhat smaller absorbing struc- 
tures) with decreasing redshift, but with a relatively small 
conversion rate of gas into stars, it might explain: (1) the 
relative flattening of the shape of the H I column density 
distribution from < z >» 2.2 to < z >« 0.8, (2) the small 
decrease in nDLA{z) per comoving volume with decreasing 
redshift, and (3) flnLA{z) showing no evidence for evolu- 
tion (or a small increase) with decreasing redshift. This 
scenario would also seem to be consistent with the mea- 
surements of DLA absorber mctallicities (§7.1) and kine- 
matics (§7.2). On the other hand, in the context of the ED 
model it is predicted that the column density distribution 
should eventually fall off as ^ N^^ at the highest column 
densities (WLFC95; Milgrom 1988). Such a fall-off would 
also be predicted if the H I gas were in randomly dis- 
tributed thin layers (Milgrom 1988). However, Figures 32 
and 33 show that the high-redshift sample does not follow 
such a trend, and that the discrepancy is even worse for 
the low-redshift sample. As mentioned above, WLFC95 
have speculated that such a fall off would take place at 
larger column densities {Nhj > 10^^ atoms cm~^) with 
increasing redshift. But with the small existing statistical 
samples of very large column density systems, it is clear 
that no such evidence for a ~ iV^^ fall-off exists. 

What is the cause of the steepening of the f{NHi,z) 
distribution for Nhi > ^0^^ atoms cmr"^ from low redshift 
to z=0? In the context of the GHC model, Khersonsky 
& Turnshek (1996) suggested that a high enough rate of 
conversion of gas into stars would naturally lead to enough 
UV radiation to destroy the more massive higher column 
density clouds, eventually causing f{NHi,z) to steepen. 
This concept is qualitatively consistent with the current 
data in that one would surmise that slsVIdla decreases by 
a factor of » 4 — 6.5 from z !v 2.0 — 0.5 to 2; = there 
must at some point be an increase in UV radiation ow- 
ing to the conversion of gas into stars. At the same time, 
the column density distribution would be expected to fall 
off as N^^ at the highest column densities in an ED 
model, and this is observed for Nhi > lO^'^ atoms cm~^ 
at z = 0. Of course, since randomly distributed gaseous 
disk galaxies were used to derive the z = column den- 
sity distribution function (Rao & Briggs 1993; Zwaan et 
al. 1999), the result at 2; = is expected. But the details 
of how the destruction of massive GHCs would finally lead 
to the column density distribution observed at z = are 
unclear. Nevertheless, the lack of a high column density 
fall-off following the ~ iV^^ functional form in the low- 
redshift DLA sample provides independent evidence that 



at < z >« 0.8 the DLA absorbers we have detected do not 
predominantly arise in sight-lines through randomly ori- 
ented gaseous disks. 

The above discussion offers some interesting possibilities 
for interpreting the evolution of the f{NHi, z) distribution 
from high redshift to z = 0. However, while it is clear that 
evolution is being observed, we cannot as yet really hope 
to use the results to provide too many details of this evolu- 
tion due to the small number statistics. Moreover, detailed 
modeling of the f{NHi,z) distribution would be compli- 
cated by the possibility that the observed f{NHi,z) dis- 
tribution arises from a population of DLA absorbers that 
has a range of properties, as opposed to sight-lines through 
randomly oriented identical absorbers. Indeed, the evi- 
dence for evolution in the shape of the f{NHi, z) distribu- 
tion suggests that the population of DLA absorbers must 
have a variety of properties. Only at z = do we observe 
clear evidence for the ~ N^^ fall-off at high column den- 
sities, consistent with the gaseous disks of spiral galaxies. 
Owing to the lack of evidence for an ^ A^/// fall-off at high 
column densities, our results on f{NHi, z) at high and low 
redshifts indicate that these samples of DLA absorbers are 
not exclusively derived from a population of randomly ori- 
ented gaseous disks, apparently ruling out the paradigm 
that DLA systems can be used to trace the evolution of 
the gaseous disk components of galaxies from high to low 
redshifts. 

6. THE LIMITATIONS OF THE RESULTS 

While we believe that the conclusions drawn thus far 

have been relatively conservative, there should be concerns 
about the reliability of some of these results and the future 
application of this technique. Some important issues are 
noted below. A clear understanding of these issues will 
require larger sample sizes and further study. 

6.1. Concerns about Small Number Statistics 

Although the results presented here represent the best 
assessment so far of DLA statistics at z < 1.65 from a 
systematic siirvey, the mimbcrs are still very small. It 
should be clear that a much larger low-redshift DLA sur- 
vey would have to be performed in order to accurately 
follow the cosmological evolution of neutral hydrogen gas 
at low redshifts. 

Our estimate of the error in the incidence of DLA in 
two redshift bins, nDLA{z ~ 0.5,1.2), as reported in §4 
was relatively straightforward. We used the Poissonian er- 
ror corresponding to the small number of detected systems 
and combined this with the error on nMgii{z) estimated 
from the SS92 Mg II survey to estimate the final error. 
However, an associated serious problem involves the de- 
termination of ^dla{z ~ 0.5, 1.2) and the corresponding 
errors. The reported floLA and errors at these redshifts 
(Figure 31) were derived using equations (9) and (10) in 
§4. However, for consistency we must also have 

^^^^ = ^^1^ / ^^^/(^^^' ^)'^^^^- 

Nhi=2x10^o 

(14) 

Thus, the value of ^Idla depends on N^f'^ . When the 
column density distributions are parameterized by power 
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laws, floLA will diverge if the power-law exponent is /3 < 2 
and N^j^ is permitted to be arbitrarily large. This is one 
reason why it has been tempting to fit the column den- 
sity distribution function with an exponential-disk model 
in the past. However, by doing so, the possible uncer- 
tainty in ^Idla{z) is minimized based on the speculation 
that we understand the model for DLA absorbing struc- 
tures. Even in a power-law model with a cutoff at N^f^ 
for /{Nffi, z), it might be reasonable to expect that N^f^ 
depends on redshift. Thus, depending on the chosen value 
of N'^f^, there could be a systematic error in the derived 
^dla{z) that is redshift dependent. As an example of 
the uncertainty, we point out that for a power law with 
index /? ss 1.7, the value of floLA increases by a factor of 
« 2.5 when N'^j^ is increased from 5 x 10^^ atoms cm^^ 
to 5 X 10^^ atoms cm~^. In DLA surveys the largest H I 
column densities are ps 5 x 10^^ atoms cm~^, so this is 
the value we adopt for the Njjf^ at all redshifts. But, 
for example, Braun (1997) has reported values as large 
as « 4 X 10^^ atoms cm^^ in local spirals. Of course, 
the actual value of N^f^ in DLA systems is very poorly 
determined because of small number statistics, and this 
problem is amplified because the highest column density 
systems are the most rare. Until sample sizes are signifi- 
cantly larger, we cannot rule out the possibility that this 
effect is causing us to underestimate ^^dla by a factor of 
as much as w 2 — 3 in a redshift dependent way. 

6.2. The Possibility of Gravitational Lensing Bias 

It has been argued that low-redshift surveys for DLA 
lines using bright QSOs may suffer a gravitational lensing 
bias in the sense that the incidence of DLA absorption 
will be systematically overestimated (Smette et al. 1997). 
This would only happen for massive lensing galaxies. For 
determination of results on QoLAiz < 1.65) this effect is 
potentially important. For example, if the two highest 
column density DLA systems contributing to the lowest 
redshift bin in Figure 31 were due to gravitational lensing 
bias, it would cause QoLAiz < 0.8) to be severely over- 
estimated. While the QSOs exhibiting low-redshift DLA 
lines observed with HST-WFPC2 do not show any evi- 
dence of being lensed (Le Brun et al. 1997), the ones that 
contributed the most to the high ri£)L^(z < 0.8) values 
have not been imaged with HST, and so the possibility of 
lensing by the very high column density low-redshift DLA 
galaxies in our survey has not been studied with the great- 
est possible sensitivity. On the other hand, our ground- 
based studies suggest that the luminous counter-parts of 
low-redshift DLA absorbers are often not highly-lmninous, 
and therefore massive, galaxies (e.g. Rao & Turnshek 
1998; §7.3) 

6.3. Reliance on Mg II Statistics 

The derived statistical results on low-redshift DLA ab- 
sorbers depend on our statistical understanding of the 
Mg II population of absorbers. At the present time the 
statistical errors on the incidence of Mg II absorption-line 
systems arc k, 47% ai z k, 0.5 and « 58% at z « 1.2. 
In the future, if we are to expand and further develop 
this technique, larger systematic surveys for QSO Mg II 
absorption- line systems (especially with Wq*^^^^> 0.5 A) 
will have to be undertaken. 



6.4. Complications Due to the Presence of Dust 

Fall & Pei (1989) originally suggested that dust in high 
column density absorbers had some influence on the statis- 
tics of DLA absorption, although at that time the obser- 
vational evidence was not very convincing. More recently, 
Pei et al. (1999) discussed the effects of dust on ^dla{z)-, 
but used the less reliable low-redshift data point of LWT95 
to show that a significant amount (»i68%) of neutral gas 
is being missed by DLA surveys. Since we derive a higher 
value for Q.dla{z) at low redshift, this effect may not be 
as severe. A comparison with their models suggests that 
we might be missing « 11% of the neutral gas mass at 
z « 0.5. 

Pettini et al. (1997a) have used depletions of Cr relative 
to Zn in DLA absorbers and found evidence for dust, but 
after taking into account the normally low metallicitics of 
DLA systems the inferred dust-to-gas ratios in the known 
population of DLA absorbers are « 1/30 of the Galaxy 
value. On the other hand, Carilli et al. (1998) have iden- 
tified some DLA absorbers by finding 21 cm absorption in 
a sample of 15 red quasars that were 1 Jy flat-spectrum 
radio sources. Based on their success, they conclude that 
optically-selected samples may miss « 1/2 of the high col- 
umn density QSO absorption-line systems due to obscura- 
tion by dust. Thus, there is reason for concern, especially 
since the effect is likely to be redshift dependent. The 
distribution of dust-to-gas ratios in high column density 
systems must be investigated in greater detail before we 
can hope to understand how dust affects our results on 
udla and VIdla- 

6.5. The HI Colurrm Density Distribution at z = 

The H I column density distribution at 2; = 0, and hence 
riDLAiz = 0) (see equation 12), was determined using ra- 
dial H I column density profiles of a complete sample of 
27 large-diameter gas-rich galaxies (Rao & Briggs 1993). 
At z = 0, spirals afford the largest cross-section to line- 
of-sight absorption at column densities Nhj > 10^^ atoms 
cm~^ (see figure 5 in Rao & Turnshek 1998) and contribute 
>95% to the local interception probability. Moreover, w 
89% of the H I mass in the local Universe resides in spiral 
galaxies later than type SO (Rao & Briggs 1993). Thus, 
the z = "DLA" statistics should be accurately described 
by the properties of local spiral galaxies. This is equivalent 
to saying that pencil-beam QSO surveys for DLA at z = 
will most likely intercept large diameter spirals. Thus, the 
statistics derived from such a survey should reproduce the 
properties of the local spiral galaxy population. 

The one caveat in using 21 cm emission measurements 
to describe column density distributions within galaxies is 
that spatial resolution is limited by the finite size of the 
radio beam. The technique measures the total H I mass 
within the radio beam to a high degree of accuracy, but 
what one really derives is an average H I column density 
within the area sampled by the beam. Thus, information 
about the highest column densities residing in any knots 
or clumps smaller than the radio beam would be lost. The 
extent to which this affects the z = statistics depends on 
the total cross-sectional area of any high-column-density 
clumps and how much smaller they are than the beam size. 
If such an effect is important, the most significant conse- 
quence might be an underestimation of f{NHi) at column 
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densities Nhi > 10^^-^ atoms cm"^ (Zwaan et al. 1999), 
but any efFcct on n]jij\(z = 0) is likely to be small since 
the highest column densities would still contribute only a 
fraction of the total cross-section over all Nhi > 2 x 10^° 
atoms cm~^. Note that the finite radio beam size does 
not affect the determination of il^LAiz — 0) since the in- 
tegrated H I mass within the beam is a well determined 
quantity. The extent to which high column densities are 
under-represented in 21 cm emission surveys can be es- 
timated by comparing the values of ilD^Aiz = 0) deter- 
mined by (1) directly measuring the H I mass of galaxies 
as done by Rao & Briggs (1993) and (2) integrating the 
/{Nhi) distribution over all column densities (see equa- 
tion (9) in WLFC95 and equations (8), (9) and (10) in 
Rao & Briggs 1993). The first is a true measure of the 
total H I mass including the contribution from the high- 
est column density clumps. If the exclusion of the highest 
column density material in the /{Nhi) distribution were 
severe, then the second estimate of the total mass would 
be significantly lower. However, this comparison was car- 
ried out in Rao & Briggs (1993) and the two values were 
found to be equal within the errors. Therefore, we con- 
clude that any high column density gas that may have 
been missed in 21 cm emission surveys is not significant, 
and that the /{Nhi) distribution at z = in Figure 33 is 
a fair representation of the local gas distribution. 

A related point involves H I self-absorption within a 
galactic disk. We note that if the inclination of a disk 
is large, the measured H I fiux in a radio beam might be 
affected by H I self-absorption within the disk. This prob- 
lem would underestimate both the mass and the mean col- 
umn density measured within the beam by approximately 
15%. We first discussed this in RTB95. While we took 
this effect into account by increasing the upper error bar 
on Q,dla{z = 0) by 15%, we did not include errors in the 
measured values of Nhi versus radius for our sample of 
spiral galaxies in the determination of /{Nhi, z = 0). 

Finally , there is the issue of small H I masses (< 
2 X 10^/ig5^ Mq) and the inability to detect them in lo- 
cal H I surveys (Zwaan et al. 1997). This is a limitation 
of our current knowledge of the distribution of H I in the 
local Universe. Observationally, the existence of absorbing 
objects with such masses (e.g. gas-rich dwarfs) that also 
contain high column density regions {Nhi > 10^^ atoms 
cm~^) cannot be ruled out. 

7. THE NATURE OF THE DLA ABSORBER POPULATION 

It is reasonable to assume that all DLA absorbers are as- 
sociated with some sort of galaxy since there are no known 
DLA systems which lack metal-line absorption. Thus, 
we can infer the presence of star formation, stellar nu- 
cleosynthesis, and stellar mass loss. In the GHC model 
for DLA systems (Khersonsky & Turnshek 1996) noted 
earlier (§5.5), it was assumed that DLA systems are asso- 
ciated with the neutral gas in GHCs that must naturally 
be associated with various types of gas-rich galaxies and 
proto-galaxies. Such galaxies would be the progenitors of 
different types of galaxies (ellipticals, spirals, irregulars, 
dwarfs, etc.) that are observed at the present epoch. This 
is consistent with the idea that the DLA systems gener- 
ally track neutral gas in, for example, both spheroidal and 
disk components of galaxies and proto-galaxies (and pos- 



sibly other gaseous structures as well) from high redshift 
to the present epoch. Alternatively, the ED model is much 
more specific in that it relates DLA absorbers to galactic 
disk systems, which furthermore are thought to be rel- 
atively luminous (WTSC86; Wolfe 1988; Lanzetta et al. 
1991; Steidel 1995). Initially there were many good indi- 
rect reasons to adopt the ED model (see Wolfe 1995 and 
references therein). However, now there are some direct re- 
sults which contradict the single DLA-H I disk paradigm. 
When evidence from DLA metallicity measurements, kine- 
matic measurements/simulations, and the direct imaging 
of galaxies associated with DLA systems is considered, the 
DLA-GHC model is seen to be a much more favorable 
framework for interpretation of the DLA systems. Here 
we summarize this evidence. The evidence supports the 
conclusion, drawn from analysis of DLA column density 
distributions at low and high redshifts (§5), that DLA sys- 
tems do not exclusively trace the gaseous disk components 
of galaxies. Instead they trace a mixture of types. 

7.1. Evidence from Measurements of the Gas-Phase 

Metallicity 

Pettini and collaborators (e.g. Pettini et al. 1997b, 
Pettini et al. 1999, Pettini et al. 2000) have measured 
metallicities in many of the high-redshift DLA systems 
and some of the known low-redshift ones. Some of their 
measurements pertain to systems which are not classi- 
cally damped, with Nhi > 2 x 10^" atoms cm^^, but 
this does not affect their broad conclusions. The most ex- 
tensive measurements are for Zn which is a particularly 
good tracer of the gas phase metallicity since it is not 
significantly depleted on to dust grains. While detailed 
interpretations of the metallicity measurements have not 
reached a consensus (e.g. see Lu et al. 1996, Pettini et al. 
1999, Prochaska & Wolfe 1999, and references therein), 
the broad results are clear. The typical column density 
weighted metallicity in the redshift interval 3 < z < 0.5 
is « 1/12 of the solar value (<[Zn/H]>Ps —1.1], with no 
obvious tendency for the metallicity to increase with de- 
creasing redshift. However, there is considerable scatter 
in the individual metallicities which is most likely caused 
by the wide range of formation histories and galaxy types 
responsible for the DLA systems. In fact, Pettini et al. 
(1999) have concluded that the known DLA systems are 
unlikely to trace the galaxy population responsible for the 
bulk of the star formation. Metallicity measurements of 
the lowest-redshift systems in our low-redshift DLA sam- 
ple will greatly clarify the evolution of DLA metallicity to 
the lowest redshifts that are currently attainable. 

7.2. Evidence from, Measurements of the Kinematics 

High-resolution studies of the metal-line profiles associ- 
ated with DLA systems provide constraints on kinematic 
models of the absorbing region. In general, a DLA ab- 
sorbing region will include one or more very high column 
density components as well as low column density compo- 
nents that give rise to detectable metal lines. Using Keck 
HIRES observations, Prochaska & Wolfe (1997; 1998) have 
shown that the kinematics of the metal-line profiles as- 
sociated with high-redshift DLA absorbers often have a 
"leading-edge" characteristic. Based on simulations, they 
have argued that the kinematic profiles of the DLA absorb- 
ing regions as a whole are more consistent with models of 
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rotating H I disks than any other single type of modcL 
More recently, the kinematics have also been shown to 
be consistent with gas in- fall due to merging (Haehnelt, 
Steinmetz, & Rauch 1998), randomly moving clouds in a 
spherical halo (McDonald & Miralda-Escude 1999), and 
multiple gas discs in a common halo (Mailer et al. 2000). 
Therefore, it seems most likely that, consistent with other 
current findings, the kinematics of DLA absorbing regions 
arise from a mix of kinematic structures, including some 
rotating gaseous disks. 

7.3. DLA Galaxy Imaging Studies 

There are some direct results from published HST- 
WFPC2 images of confirmed {Nhi > 2 x 10^° cm'^) low- 
redshift DLA absorbers which suggest that the morpho- 
logical types of DLA galaxies are indeed mixed. Seven 
DLA absorbers in the redshift interval 0.395 < z < 1.010 
have so far been imaged with _ffS'T-WFPC2. In one case 
(the z = 0.656 system toward 3C 336 (1622-^239)) there 
is no optical evidence for a galaxy despite very deep ob- 
servations (Steidel et al. 1997). In the other six cases the 
galaxies are likely to be spirals (the z = 0.437 system to- 
ward 3C 196 (0809-F483) and the z = 0.633 system toward 
Q1209-I-107), two amorphous, LSB galaxies (the z — 0.395 
system toward PKS 1229-021 and the z = 0.692 system 
toward 3C 286 (1328-1-307)), and two compact objects (the 
z = 1.010 system toward EX 0302-223 and the z = 0.860 
system toward PKS 0454+039) (Lc Brun et al. 1997). One 
other higher-redshift DLA absorber imaged by Le Brun et 
al. (1997) with i7S'T-WFPC2 was also found to be a com- 
pact object (the z = 1.776 systems toward MC 1331-1-170). 

In addition, ground-based observations of three of the 
lowest redshift DLA absorbers in our low-redshift sample 
have been published. Initial WIYN imaging indicates that 
one is a dwarf while the other is a dwarf and/or an object 
with very low surface brightness (the z — 0.091 and 0.221 
absorbers toward 01 363 (0738+313)) (Rao & Turnshek 
1998). In the third luminous spiral lies near the 

sight-line but several dwarfs are also present with even 
smaller impact parameters (the z = 0.313 system toward 
PKS 1127-14, Lane et al. 1998a). 

Thus, the DLA galaxies associated with the DLA sys- 
tems in our low-redshift sample are not dominated by 
luminous gas-rich spirals. Instead, the imaging results 
demonstrate that the morphological types of DLA galax- 
ies arc mixed and span a range in luminosities and surface 
brightnesses. Ongoing imaging observations of other DLA 
absorbers in our low-redshift sample support this conclu- 
sion. 

7.4. Discussion 

Several independent lines of evidence now lead us to con- 
clude that the DLA galaxy population represents a range 
of morphological types and luminosities. It might then 
follow that the progenitors of spiral galaxies are not the 
main repositories of neutral gas at redshifts z > 0. How- 
ever, this is at odds with results near redshift z = 0, where 
the bulk of the high-column-density neutral gas is found 
to reside in luminous gas-rich spirals. Can these two ap- 
parently disparate observational results be reconciled? 

When considered together, the following three results, 
(1) that riDLAiz) and i^oLAiz) show little evolution be- 
tween redshifts 2: w 4 and z ~ 0.5, (2) that the metallici- 



ties of the DLA absorbers has remained roughly constant 
over this redshift interval, and (3) that many of the DLA 
galaxies at low redshift that have been identified through 
imaging are LSB or dwarf galaxies, might lead us to con- 
clude that DLA absorbers are a slowly evolving popula- 
tion. In fact, DLA galaxies probably do not trace star- 
forming galaxies (see also Fynbo et al. 1999), and per- 
haps DLA absorbing regions are in a quiescent phase of 
their evolution when they are not undergoing an episode 
of star formation. Thus, we might be seeing little evo- 
lution in the properties of DLA systems by virtue of an 
observational bias that selects similar types of absorbers 
at all redshifts, i.e., regions in intervening galaxies with 
high-column-density, low-metallicity (and hence, low-dust- 
content) gas. If the models of Pei et al. (1999) are correct, 
then galaxies that are not tracked by DLA absorption lines 
make up w 10% of the neutral gas at 2; ~ 0.5 (using our 
new low-redshift result) and « 60% at z « 2.5 where the 
effect of obscuration by dust is most severe (see their fig- 
ure 7). On the other hand, local H I 21 cm cnnission-linc 
surveys do not suffer from this bias. Thus, we are left with 
some unanswered questions: (1) Are a majority of spiral- 
galaxy progenitors missing from DLA absorption-line sur- 
veys and do they contain the missing gas fraction seen in 
the models of Pei et al. (1999)? (2) What fraction of the 
neutral gas mass is contained in star-forming galaxies at 
high redshift? (3) What are the present-day whereabouts 
of the neutral gas observed in high- and low-redshift DLA 
galaxies? 

A self-consistent explanation of the observational data 
that has been obtained from QSO absorption lines and 
imaging studies of DLA galaxies, UV luminous galaxies, 
faint blue dwarfs, etc. is still far from complete. Re- 
call that the redshift interval < 2; < 0.5 includes the 
most recent 45% of the age of the Universe for qq = 0.5 
(33% for go = 0). Thus, it is conceivable that significant 
changes in the properties of galaxies have taken place over 
this time interval. Efforts on increasing the sample size 
of low-redshift DLA absorbers should be made. This will 
allow us to study the evolution in nDLA{z), ^dla{z), and 
/{Nhi, z) in greater detail, and help assess the effects of 
biases such as dust obscuration. 

8. SUMMARY AND CONCLUSIONS 

We have used a novel technique, based on Mg II pre- 
selection, to conduct a low-redshift {z < 1.65) survey for 
DLA absorption systems in QSO spectra. The survey was 
conducted in the UV with HST-FOS. Our low-redshift re- 
sults, when combined with results from DLA studies at 
higher redshift, afford us the opportunity to study the dis- 
tribution of neutral gas in DLA systems back to a time 
when the Universe was less than 10% of its present age. 
Aside from DLA studies, no other methods have revealed 
such large amoimts of significantly redshifted neutral gas. 
Thus, the study of DLA absorbers is currently the most 
effective method for investigating the bulk of the known 
neutral gas mass of the Universe. Our findings and con- 
clusions are as follows: 

(1) In our survey we uncovered 12 DLA lines in a sam- 
ple of 87 Mg II systems with Mg II rest equivalent width 
I^A2796> 0.3 A. Two more DLA systems with no avail- 
able Mg II absorption-line information were discovered 
serendipitously. 
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(2) The DLA absorbers are drawn almost exclusively 
from the population of Mg II absorbers which have 
I^A2796> Q_g ^_ Approximately 50% of the systems with 
Mg II A2796 >0.5 A and Fe II A2600 >0.5 A have DLA 
absorption lines which meet the classical definition used 
in high-redshift surveys. Nhi >2x 10^" atoms cm~^. We 
know of only one rare case of a low-redshift DLA absorber 
that does not have Mg II A2796 >0.5 A and Fe II A2600 
>0.5 A. 

(3) The survey method resulted in a determination of 
the fraction of Mg II systems with DLA, which in turn 
allowed us to determine the incidence of low-redshift DLA 
absorption using earlier results on the incidence of Mg II 
absorption systems. Neutral hydrogen column densities 
were determined by fitting Voigt damping profiles to the 
Lya lines in the HST-FOS UV survey spectra. This in- 
formation was then used to deduce the low-redshift cos- 
mological neutral gas mass density of DLA absorbers and 
their H I column density distribution. 

(4) In making these determinations, possible biases were 
carefully eliminated. For example, known 21 cm absorbers 
were excluded from the initial low-redshift sample and the 
low-redshift sample (this study, < z >« 0.8) and high- 
redshift sample (from the literature, < z >« 2.2) were 
corrected for Malmquist bias before they were compared. 
A comparison with inferred results on DLA at 2; = was 
also made; the z = results were derived from 21 cm 
emission observations of nearby gas-rich spirals, since spi- 
rals are the dominant known reservoirs of large column 
density neutral gas in the local Universe. 

(5) The findings were discussed in the context of a A = 
cosmology with qq = and = 0.5. Assuming that the 
z = results arc not missing any significant fraction of 
the local neutral gas, the incidence of DLA systems per 
unit redshift, njjLA, is found to decrease with decreasing 
redshift by an amount significantly larger than what is ex- 
pected in a no-evolution universe. Specifically, we found 
that the high-redshift incidence of DLA systems is w 6 
times larger than that found at z = in a go = no- 
evolution universe (» 10 times larger if go = 0.5). How- 
ever, if the z = point is removed from the analysis, there 
is only mild evidence for a real decrease in incidence per 
comoving volume down to z ^ 0.5. Thus, it is possible 
that most of the real evolution occurs between z « 0.5 
and z = 0. 

(6) The DLA cosmological neutral gas mass density also 
shows strong evidence for evolution when compared with 
the z = results. If qo = 0, the comoving cosmological 
mass density of high-redshift DLA absorbers is w 4 times 
larger than what is inferred from local 21 cm observations 
(« 6.5 times larger if qq = 0.5). However, if the z = re- 
sults are ignored, the cosmological mass density of neutral 
gas in low-redshift DLA absorbers is observed to be com- 
parable to that observed at high redshift (and, in fact, for- 
mally somewhat higher). Thus, we cannot point to a trend 
in the DLA data themselves which indicates that ^Idla at 
low redshift is approaching the value derived for the cos- 
mological mass density of neutral gas found in present-day 
spirals. Again, it is possible that most of the real evolution 
occurs between z « 0.5 and z = 0. 

(7) The H I column density distribution of the low- 
redshift DLA sample is found to differ in comparison to the 



high-redshift DLA sample and that inferred for local spi- 
rals. The low-redshift DLA absorbers exhibit a larger frac- 
tion of very high column density systems in comparison to 
determinations at both high redshift and locally. The ob- 
served trend suggests that the column density distribution 
starts out relatively steep at high-redshift (/ ~ NJj]'^), it 
then flattens at the redshifts covered in our low-redshift 
survey (/ ~ N^]'^), and then locally it starts off relatively 
flat (/ ~ N^y^) but then steepens at Nhi > 10^"^ atoms 
cm^^ (/ ~ N^f'^), becoming even steeper than observed 
at high redshift. We emphasize that the column density 
distributions of DLA absorbers at high and low redshifts 
are not observed to fall off as / ^ N^^ with increasing 
column density. An / ~ N^^ fall-off is theoretically pre- 
dicted for disk-like systems and this is, in fact, what is 
observed locally in spiral samples. 

(8) While our understanding is that it is unlikely that 
the z = studies have missed any signiflcant component 
of the neutral gas, we also understand that existing H I 
21 cm emission surveys are only very sensitive to neutral 
gas masses Mhi > 2 x lO'^h^^ Mq. Thus, if the z = 
results are missing a significant component of the neutral 
gas mass, it almost certainly lies in this regime. How- 
ever unlikely, if this should turn out to be the case, there 
would then be evidence for only mild intrinsic evolution in 
the DLA population over the redshift interval 3 > z > 0. 

(9) The decrease in the overall incidence of DLA ab- 
sorbers with decreasing redshift, accompanied by the in- 
crease in the relative number of high column density sys- 
tems with decreasing redshift, explains the lack of evidence 
for evolution in the amount of neutral gas locked up in the 
DLA absorbers from high redshift to z « 0.5. This lack 
of evidence for significant conversion of gas into stars (or 
some other phase), suggests that the gas contained in DLA 
systems has not experienced much star formation, and is 
consistent with the lack of evolution seen in their metallic- 
ities. Moreover, recent results from imaging studies show 
that DLA galaxies comprise a mix of morphological types 
with a large fraction of them being low-surface-brightness 
or dwarf galaxies. Thus, DLA absorption lines seem to 
be tracking a slowly evolving population of objects. In- 
deed, there might be little evolution in the properties of 
DLA systems for the simple reason that observational bi- 
ases select similar types of absorbers at all redshifts, i.e., 
regions in intervening galaxies with high-column-density, 
low-metallicity (and hence, low-dust-content) gas. 

(10) DLA systems trace the bulk of the observable neu- 
tral gas content of the Universe at redshifts z > 0, while at 
the present epoch, the ncnitral gas mass density and cross- 
section are dominated by luminous spirals. In order to un- 
derstand the evolution of neutral gas in the Universe, we 
must address the question of the present-day whereabouts 
of the neutral gas mass observed in high- and low-redshift 
DLA systems. Can models predict the decline of a factor 
of ~ 5 in the neutral gas mass density from z « 0.5 to 
z = especially if DLA galaxies are a non-star-forming 
and slowly-evolving population? 

(11) Finally, in order to improve our understanding of 
the neutral gas content of the Universe to the point where 
we can accurately follow its evolution to z = 0, we will 
have to collect and analyze large samples of low-redshift 
DLA systems. In addition to improving the statistical un- 
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certainties, this would permit a better determination of 
the shape of the H I cohimn density distribution at large 
column densities, effectively allowing a determination of 
the maximum H I column density which might be a func- 
tion of redshift. Understanding this turn-down in column 
density is important for the determination of Qdla- A 
large optical ground-based survey for large rest equivalent 
width QSO Mg II absorption lines would be a reasonable 
initial step. Additionally, however, in order to fully take 
advantage of the information that could be provided by a 
large DLA survey, we would have to do more than sim- 
ply increase the sample size. We would also have to de- 
sign observational programs that clarify the uncertainties 
or systematic effects caused by potential dust obscuration 
and possible gravitational lensing bias. The sensitivity of 
the H I 21 cm emission surveys should also be pushed to 
lower limits, to the point where we can definitively assess 
the possibility of very high column density DLA absorp- 
tion arising in very low mass local objects {Mhi < 10^ 
Mo). 
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